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Abstract 

SUSYGEN is a Monte Carlo program designed for computing distributions and generating events 
for MSSM sparticle production in e + e~ collisions. The Supersymmetric (SUSY) mass spectrum 
may either be supplied by the user, or can alternatively be calculated in two different models of 
SUSY Breaking: gravity mediated supersymmetry breaking (SUGRA), and gauge mediated super- 
symmetry breaking (GMSB). The program incorporates the most important production processes 
and decay modes, including the full set of R-parity violating decays, and the decays to the grav- 
itino in GMSB models. Initial state radiation corrections take into account pt/pl effects in the 
Structure Function formalism, and an optimised hadronisation interface to JETSET 7.4 including 
final state radiation is also provided. 
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NEW VERSION SUMMARY 



Title of new version: SUSYGEN 2.2 
Catalogue number: 

Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland (see appli- 
cation form in this issue) 

Reference to original program: SUSYGEN 1.0; S. Katsanevas and S. Melanchrinos, "SUSYGEN'', in 
the proceedings of the Workshop Physics at LEP 2, eds. G. Altarelli, T. Sjostrand, and F. Zwirner, 
CERN 96-01, page 328, Geneva, 19 February 1996. 

Authors of original program: Stavros Katsanevas, Stavros Melachroinos 
The new version supersedes the original program 
Licensing provisions: none 

Computer for which the new version is designed: 
DEC ALPHA, HP 9000/700, IBM AIX series; 



WWW site: http: / /lyohp5.in2p3.fr/delphi/katsan/susygen.html 

Operating system under which the new version has been tested: UNIX, VMS 

Programming language used in the new version: FORTRAN 77 

Memory required to execute with typical data: ~ 250 - 450Kb 

No. of bits in a word: 32 

No. of processors used: 1 

The code has not been vectorised 

Subprograms used: CERNLIB [i], JETSET [ii], RANLUX [Hi], PHOTOS [iv] 
No. of lines in distributed program, including test data, etc. : ~ 18000 
Correspondence to: Stavros. Katsanevas@cern.ch, Peter. Morawitz@cern.ch 

Keywords: e + e~ collisions, LEP, Super symmetry, SUSY, MSSM, R-parity violation, Supergravity, 
Gauge mediated Supersymmetry breaking 

Nature of physical problem 

The complexity and the interdependence of the MSSM signals demand a fast and precise simula- 
tion program which treats the production and decay of all possible sparticles in different models of 
supersymmetry breaking, and in models in which R-parity is either conserved or violated. 

Method of solution 

SUSYGEN 2.2 is a Monte Carlo generator for the tree-level production and decay of all MSSM sparti- 
cles at e + e~ colliders. It includes pair-production of gauginos, sfermions and squarks, the production 
of a gravitino plus a neutralino in Gauge-mediated SUSY breaking models, the single production of 
gauginos through s-channel resonance in R-parity violating models, and the production of Higgses. 
The program implements all important decay modes of sparticles relevant to LEP energies, including 
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cascade decays to lighter SUSY states, radiative decays of gauginos to a photon, decays to the Higgses, 
decays to gravitinos, and R-parity violating decays to Standard Model particles. The precise matrix 
elements are used for the production and for the decay of sparticles, though helicity correlations are 
not taken into account. The mass spectrum of the SUSY particles at the electroweak scale is either 
calculated assuming a common mass scale at the GUT scale (Supergravity models or Gauge mediated 
SUSY breaking models), or the mass spectrum can be optionally provided by the user. The LSP can 
either be the neutralino, the sneutrino or the gravitino in R-parity conserving models, or there can be 
no LSP in the sense that all sparticles can decay to standard particles through R-Parity violation. The 
initial state radiative corrections take account of Pt/pl effects in the Structure Function formalism. 
QED final state radiation is implemented using the PHOTOS library. An optimised hadronisation 
interface to JETSET 7.4 is provided, which also takes into account the lifetime of sparticles. 

Reasons for the new version 

SUSYGEN has been extensively used by the four LEP collaborations to simulate the expected signals 
on the way to the higher LEP energies. The scope of version 1.0 of the program was limited in the 
following way: only neutralino LSPs were implemented; gauge-mediated SUSY breaking scenarios, 
i.e. the production and decay of gravitinos was not implemented; no R-parity violation option was 
available; QED final state radiation was not included, and none of the Higgs production mechanisms 
were implemented. 

Restrictions on the complexity of the problem 

The factorised treatment of the production and decay of sparticles through the appropriate matrix ele- 
ments is more accurate than the equivalent treatment of generators like e.g ISAJET M and SPYTHIA 
M, where the sparticle decays are described by phase space. Nevertheless it is less accurate than pro- 
grams which take helicity correlations into account, as for example is done in SUSY-GRACE (SUSY23) 
||, or DFGT j|. Squark hadronisation before decay is not implemented, and the hadronisation of 
squark decays via the R-parity violating coupling A" is only simulated using the independent frag- 
mentation scheme. The decays to gluinos are at present not implemented. And three-body decays of 
sfermions to lighter sfermions are not included. 

Typical running time: 40 events/sec on a DEC ALPHA. 

Unusual features of the program: none 
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LONG WRITE-UP 



1 Introduction 

In recent years supersymmetry (SUSY) has been extensively used to chart the map of possible physics 
beyond the standard model (SM). This symmetry predicts the existence of additional particles which 
differ from their standard model partners by half a unit of spin. The simplest model available is 
the Minimal Supersymmetric Standard Model (MSSM) 0, ||] which contains the minimal number 
of new particles and interactions that are consistent with the SM gauge group. In the MSSM all 
SM fermions have scalar SUSY partners (see Table |l]): the sleptons, sneutrinos and squarks. The 
SUSY equivalent of the gauge and Higgs bosons are the charginos and neutralinos (the gauginos), 
which are the mass eigenstates of the (W + ,H + ) and (7, Z, H®, H®) fields, respectively. Within the 
supersymmetric framework many theoretical questions of grand unified theories (GUTs) such as the 
hierarchy problem and the unification of couplings pf] may be successfully addressed. 
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W^,!!^ mix to form 2 chargino mass eigenstates XxiXi 
7, Z, H°i t 2 mix to form 4 neutralino mass eigenstates Xi^X^X^Xi 
tL,tR (and b, f) mix to form the mass eigenstates £1,^2 



Table 1: Particle-Sparticle correspondence in the MSSM 



The masses and couplings of the new supersymmetric states are related by the symmetry to those 
of the SM states. If supersymmetry were an exact symmetry, the SUSY partners would be degenerate 
in mass with their SM partners. Since no such states have been observed to date, supersymmetry 
must be broken at some higher energy scale. Two supersymmetric breaking mechanisms have been 
studied in detail in the literature: gravity mediated supersymmetry breaking (SUGRA) || |i| and 
gauge mediated supersymmetry breaking (GMSB) p], 11]. 



The first model (SUGRA) is inspired by supergravity theories, where supersymmetry is broken at 
some very high scale, a "hidden sector", close to the Grand Unification scale and is communicated 
to the visible sector through gravitational interactions. If one assumes universal soft breaking terms 
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at the GUT scale, the masses of the SUSY particles and their couplings can be calculated at the 
electroweak scale through the evolution of the renormalisation group equations (RGEs). The model 
thus predicts the entire SUSY spectrum from a few parameters. 

Recently the second model (GMSB), where the supersymmetry breaking occurs at a relatively low 
scale (a few hundred TeV), has received considerable attention in the literature. Here the messengers 
of the supersymmetric breaking are the gauge bosons. This model also predicts mass relationships 
for the gauginos and the scalars, but the main difference to SUGRA models is that the G goldstinof] 
- the longitudinal component of the gravitino - can be very light, and would thus be the Lightest 
Supersymmetric Particle (LSP). 

Apart from the issue of supersymmetry breaking, a further complication arises when the SM is 
extended to incorporate SUSY. The most general interactions of the SM and SUSY particles invariant 
under the SU(3) C x SU(2)l x U(l)y gauge symmetry are those of the MSSM plus the additional 
superpotential terms Jl2[ 

Wij^ = XijkLiLjE k + X' i j k L i QjD k + \"j k UiDjD k . (1) 

Here L (Q) are the lepton (quark) doublet superfields, and D, U (E) are the down-like and up-like 
quark (lepton) singlet superfields, respectively; A, A', A" are Yukawa couplings, and k = 1,2,3 are 
generation indices. The first two terms in Eq. (|l]) violate lepton number, while the last term violates 
baryon number. The simultaneous presence of the last two terms leads to rapid proton decay, and the 
solution of this problem in the MSSM is to exclude all terms in Eq. (||) by imposing the conservation 



of a discrete multiplicative quantum number, R-parity [L4j], defined as 



R p = { _lfB+L+2S (2) 

Here B denotes the baryon number, L the lepton number and S the spin of the particle state. However, 
this solution is not unique, and a number of models |0| predict only a subset of the terms in ([[]), 
thus protecting the proton from decay. These alternative solutions are denoted "R-parity violation". 
Theoretically there is no clear preference for models which impose R-parity conservation or violation 
[fl4| . Experimentally the two cases correspond to vastly different signatures, as will be discussed below. 

If R-parity is conserved one may expect two different classes of candidates for the LSP, depending 
on whether the messengers of the supersymmetry breaking are gravitational or the gauge bosons: 
the neutralino (or the sneutrino), or the gravitino. In either case the LSP is stable, a consequence 
of R-parity conservation. If the neutralino (or the sneutrino) is the LSP one expects events with 
large missing energy (E ) due to the neutralino escaping detection. In the second case where the 
gravitino is the LSP, the possible signatures depend on the decay rate of the heavier SUSY particles 
to the gravitino and on the nature of the NLSP (Next Lowest Supersymmetric Particle). For example, 
neutralino NLSPs decay to an undetectable gravitino and a photon, while sfermion NLSPs decay 
to the corresponding fermion and the gravitino. Prompt decays to the gravitino would produce the 
classical missing energy signatures, possibly with additional photons in the event; slow decays inside 
the detector volume may produce tracks with kinks or displaced vertices; and decays outside the 
detector may also produce heavy stable charged particle signatures if the NLSP is a charged particle. 

R-parity violation has two mayor consequences for collider phenomenology. Firstly, the LSP is not 
stable and can decay to SM particles. And secondly, sparticles may be produced singly. This opens 
up a a whole new area of different signatures, which in the simplest case may consist of acoplanar 
lepton pairs or acoplanar jets, or in the more complex case of multi-lepton and/or multi-jet final states, 
possibly with some missing energy. 

1 Throughout this paper we will refer to the G as the gravitino. 
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SUSYGEN 2.2 is a Monte Carlo generator capable of simulating the production and decay of all 
possible sparticles in the different scenarios described above, i.e. in models of different supersymmetry 
breaking (SUGRA or GMSB), and in models in which R-parity is either conserved or violated. SUSY- 
GEN is dedicated to e + e~ colliders^ where it has been extensively used (and debugged!) by the four 
LEP collaborations. The program is based, for the most part, on the formalisms of jl?]] for the MSSM 



sparticle production and decays, on [18, 19] for SUGRA, on |2(J for GMSB, and on [22, 23] for the 
ftp production and decays. It uses the precise matrix elements for the production and for the decay 
of sparticles, although helicity effects are not taken into account. SUSYGEN incorporates initial and 
final state radiative corrections, and includes a hadronisation interface to JETSET 7.4 [24|. 

The outline of this paper is as follows: after reviewing and defining the MSSM framework in 
Section E3, the physics content of the generator is described in Section ||. The structure of the program 
is outlined in Section and instructions on how to obtain, compile and run the program are given in 
Section |5| Conclusions and future plans are given in Section ||. The Appendix lists the particle and 
production code conventions, the gaugino mass matrix basis, the formulas used for the production and 
the decay of sparticles, common blocks and the ntuple variables. Finally a test run input and output 
listing is given. 



2 The other three "general" SUSY programs are ISAJET §, SPYTHIA g] and SUSY-GRACE (SUSY23) ||. They 
have different areas of speciality: ISAJET and SPYTHIA treat both hadronic and e + e~ colliders, but they treat the 
sparticle decays through phase-space and they do not contain the full set of R p decays, while SUSY-GRACE includes 
helicity correlations, but does not incorporate all SUSY processes. It has no ft p option yet. See also |l6| and references 
therein for other generators or programs related to SUSY. 
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2 The MSSM framework 



The MSSM is defined as the Supersymmetric extension of the SM which contains the minimal number 
of new particles and interactions that are consistent with the SM gauge group. The gauge couplings 
of the SUSY particles are equal to the SM gauge couplings. Experimental evidence shows that SUSY 
must be broken, and the masses of the superpartners are therefore much heavier than their SM 
counterparts. The exact nature of the Supersymmetry breaking mechanism is unknown, but it can be 
parameterised by the following soft supersymmetric breaking terms p5[ : 

• Gaugino masses: MjAjAj, where i = 1,2,3 for U(l), SU(2), SU(3) and A» denotes the gaugino 
partner of the corresponding gauge field. 

• Higgsino masses: [ihxh® and fihfh^, where /x is the higgsino mass mixing parameter. 

• Scalar masses: e.g. (^l^l + ^ip£) + m \^R^R + m \ ^r^R, plus similiar expressions for the 
first two generations and the sleptons and sneutrinos. 

• 'A' Left-Right trilinear mixing terms: e.g. AtXtithH® — bj^H^)^, plus similiar expressions for 
At and A T . 

In the "unconstrained MSSM" the above mass terms are considered to be free parameters of the 
theory (Section |2.1|) . The large number of unknowns can be substantially reduced in "constrained 
versions of the MSSM" if a particular SUSY breaking mechanism is employed. Two such mechanisms 



are implemented in SUSYGEN, SUGRA and GMSB, and are discussed in Section |2.2| . One of the most 
important aspects for collider phenomenology is the nature of the LSP, which determines the dominant 
decay modes of SUSY particles and hence their experimental signatures. This will be discussed in 



Section 2.3 



2.1 The unconstrained MSSM 

The unconstrained MSSM corresponds in SUSYGEN to the following choices of the input card 



"MODES" (see also Section |5j) : MODES=2, MODES=4 or MODES=5. 



2.1.1 Gaugino Masses 

There are four spin 1/2 SUSY partners to the colorless neutral gauge fields and higgs bosons, the 



neutral "gauginos". Their interaction eigenstates may be written as [17] 



$ = (-«7, -iZ, cos PH% - sin /3iJ°, sin pH$ + cos /3H$) , (3) 

where 7, Z, H®, denote the the two component spinor fields of the photino, zino and the two neutral 
higgsinos, respectively; sq and cq are the weak angle sine and cosine; and tan/3 = V1/V2, the ratio of 
the vacuum expectation values of the neutral Higgs fields. The mass matrix of the neutral gauginos 
is given by 

( M 1 c 2 e + M 2 s 2 (M 2 -M 1 )c e s \ 

(M 2 -Mi)c 9 s M 2 c 2 g + M lS 2 e m z 

mz H sin 2/3 —/i cos 2/3 

\ -/icos2/3 -/isin2^ J 

The physical mass eigenstates, the neutralinos (x^)> are defined by 



(4) 



Xk = N kl ip? (5) 
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and the neutralino mass eigenvalues (m x o ) can be computed by diagonalising the mass matrix Eq. (||) 
according to 

m x o5 k i = N km N ln M^ n . (6) 
There are two SUSY partners to the W gauge boson and the charged Higgs bosons. Their mass 



eigenstates are the charginos (xt > X k )> which are admixtures of their interaction eigenstates [17] 

= (~iW + ,H+) 
If = (~iW-,H{). (7) 

The chargino mass matrix is given by 

M c _ ( M 2 \/2m w sin \ . g . 

I \f2mw cos (3 /j, I 

and the chargino mass eigenstates are given by 

xt = W 

X^ = U kl iP~ (9) 

The chargino mass eigenvalues (m +) can be computed by diagonalising the mass matrix Eq. (R) 
according to 

m x +S kl = U* km V ln M^ n . (10) 

The diagonalisation of the above two mass matrices gives four neutralino (x?, X% X% X4) an d two 
chargino (Xx 1X2) physical states. The masses and the couplings of these states are determined by the 
eigenvalues and the eigenvectors of Eq. (|0|,[lO|). They depend on the three parameters (M 2 , tan/3, /i) 
in the case of charginos, and on the four parameters (Mi, M 2 , tan/9, /i) in the case of neutralinos. 
Appendix |B| gives more details on gaugino mixing and the different neutralino bases used inside 
SUSYGEN. 

In order to simplify the parametric dependence a relationship between the gaugino masses is 
commonly assumed. In the simplest possible picture the gaugino masses Mi,M 2 and M3 are assumed 
to be equal at the grand unification scale, an assumption which is natural in the context of gravity 
induced supersymmetry breaking. The one-loop renormalisation group equations then allow one to 
calculate the evolution of the three couplings to the EW scale. This determines the "GUT relations" 
for the gaugino mass terms [25]: 

Mi = ^tan 2 ^M2 (11) 
M 3 = -^sm 2 9 w M 2 (12) 

OEM 

Numerically this relationship (at the EW scale) is approximately given by M3 : M 2 : M\ ~ 7 : 2 : 1. 
The above relationship obtained in the context of SUGRA models is also valid for GMSB models. 
However, the existence of string models or other models with non-universal masses at the GUT 



scale [26] induces some scepticism on the general validity of Eq. ( |lT| , |l^ ) . By default SUSYGEN uses 
Eq. (11), but the ratio M2 : Ml may be optionally changed by the user. 

The GUT relations ensure that the gluino (where m.g = M3, neglecting QCD corrections) is 
normally much heavier than the SU(2) and U(l) gauginos. Their production is therefore typically 
inaccessible at LEP. However, in some models loop corrections to the gluino masses predict light 
gluinos [27|. And in GMSB models with two effective SUSY breaking scales, one for a colour-singlet 



and one for a colour non-singlet sector, gluinos may also be light, possibly even the LSP 28]. Light 
gluino scenarios are at present not implemented in SUSYGEN. 
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2.1.2 The Gravitino Mass 

In GMSB models spontaneous SUSY breaking leads to a goldstone fermion, the goldstino G with very 
low mass. In local SUSY G is the longitudinal component of the gravitino. The G mass is determined 
by F, the characteristic scale of SUSY breaking, and is given by [25]: 



= = 2.5 x ( ff— feV (13) 

For y/F ~ 10 8 TeV - the characteristic scale of string-motivated supergravity models - the G is fairly 
massive and has no consequences for LEP phenomenology. In gauge mediated models \~F ~ 100—2000 
TeV, and the G is the LSP. Note that the scale \/F also determines the coupling strength of the 
gravitino and thus enters in the neutralino-gravitino production cross section, and in the decay rate 
of the lifetime of the NLSP, which can only decay to the gravitino if R-parity is conserved. 

2.1.3 Third Generation Sfermion Masses and Mixing Angles 

The SUSY partners of the SM fermions are the sleptons Il,Ir, the sneutrinos v and the squarks qi,QR- 
The charged leptons and the quarks have two SUSY partners, one for each lepton/quark chirality. The 
third generation sfermions (t, 6, f ) are the most likely candidates for the lightest scalar quark or lepton 
states because of the potential for large mixing angles between the left and right handed states, and 
because of the large third generation Yukawa couplings. This may be seen from the stop and sbottom 
mass matrices fl2l 




m \ L + m t + m i) rntXt 

mtXt m \ R ~^ m t m< D 

m\ + m? + m 2 n rn^Xf, 

OL 

mbXb m~ + ml + m 2 D 




(14) 



from which the mass eigenstates of the stop (£1,^2) and the sbottom (61,62) can be obtained by 
diagonalisation. Here 

X t = A t -ficot(3 (15) 

X b = Ab-jutan/? (16) 

m 2 D = (I 3 - Q S in 2 9 W )\ cos 2 (3\m 2 z (17) 

where mo are the so-called "D-terms", I 3 is the third component of isospin, and Q is the electric 
charge. Note that for stops and sbottoms the off-diagonal terms in Eq. ( |l4| ) are sizeable and lead 
to appreciable mixing between the left-right states. This is not the case for the first two squark 
generations, where the analogous m q X q terms are small since m q ~ 0. Depending on tan (3 either the 
stop or the sbottom can be the lightest squark state. A similar argument applies for the stau, which 
can be lighter than the selectron and the smuon due to mixing. 

SUSYGEN implements mixing between the left-right third generation sfermions, and the user 
can control the amount of mixing through the parameters At, Ab, A T . In the unconstrained MSSM 
(MODES=2,4) the user may alternatively provide the explicit mixing angles 9^, 9^, 9f which parame- 
terise the mass-eigenstates in terms of the left-right states: 

/1 = h cos 0j + f R sin 9 j 

h = -f L sm9 I + f RC o S 9 I (18) 
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2.1.4 The Higgs Mass 

The tree level Higgs masses only depend on two parameters: tua, the mass of the CP-odd Higgs and 
tan j3. The radiative corrections |3(J modify these predictions and introduce a dependence on the rest 
of the supersymmetric mass spectrum. The RGE improved formulae [19] implemented in SUSYGEN 



take into account most of these corrections. The Higgs masses show a strong dependence on the mixing 
of the stop Xt, increasing the theoretical upper limit on the Higgs mass for larger values of Xf. 

2.2 The constrained MSSM 

The constrained MSSM corresponds in SUSYGEN to the following choices of the input card "MODES" 
(see also Section MODES=l or MODES=3. 

2.2.1 Mass Spectrum in SUGRA models 

In Supergravity inspired models the soft SUSY breaking parameters are assumed to be universal at 
the GUT scale, reducing the number of parameters to: 

• mo, the common mass of the sfermions at the GUT scale. 

• M%, the SU(2) gaugino mass. 

• /x, the mass mixing parameter of the Higgs doublets. 

• tan (3, the ratio of the vacuum expectation values of the two Higgs doublets. 

• At, Ah, A T , the trilinear couplings in the Higgs sector. 

All above parameters are defined at the EW scale except for tuq, which is defined at the GUT scale^. 
The sfermion masses tuq are evolved from the GUT scale to the EW scale according to the formulae 



given in Appendix B of reference [18]. Mixing of the third generation sfermions is taken into account 
according to Eq. ( 14 ) through the parameters At , A^ and A T . For certain sets of parameters the 
sfermion masses at the EW scale may turn out negative, and in this case a warning is issued and 
the program skips this particular parameter point. The Gaugino and the Higgs mass spectra are 



calculated as discussed in Sections |2.1.l]]2,l~4 



For the first two generations, and neglecting the D-terms Eq. (|17|), the renormalisation of the scalar 
masses at e.g. tan/? = 1 give 

(m| - ml) : (m| - m§) : (mf - mg) : M 2 = 0.22 : 0.75 : 6.40 : 1 (19) 

displaying a clear mass hierarchy, with the right sleptons (mjj) being the lightest states, followed by 
a mass degenerate left slepton doublet (ra^), and the first and second generation squarks at higher 
masses. 

Radiative corrections to the soft SUSY breaking terms can induce electroweak symmetry breaking, 
resulting in further constrains on the parameters of SUGRA models. This is at present not imple- 



mented within SUSYGEN, although private versions of SUSYGEN interfaced to the code of [31] exist, 
and may be obtained from the authors. 

2.2.2 Mass Spectrum in GMSB models 

In GMSB models the scalar mass relationships are p5|]: 



mi : ml : m 2 R : Ml = 11.6 : 2.5 : 1.1 : JN 5il0 (20) 



3 Note that for example the ISAJET program Q defines the above parameters at the GUT scale. 
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Here N^xo is a number characterising the number of representations in the messenger sector. Note 
that the squarks are in principle very heavy in GMSB models. The gravitino is the LSP, and its mass 
is given by Eq. (|i~3|). For N^xo = 1 the neutralino is the NLSP, while for N$xo — 2 the Ir is the 
NLSP. The Gaugino and the Higgs mass spectra are calculated in the same way as in the SUGRA 
model. The following parameters determine the entire mass spectrum of the model: -/V510, fJ,, M2, 
tan/3, A t ,A b ,A T . 



2.3 The Nature of the LSP and the Sparticle Decay Modes 

SM particles have R p = +1 while their SUSY partners have R p = —1, which can easily be seen from 
Eq. (0) . The multiplicative conservation of R-parity only allows vertices with an even number of SUSY 
particles, and the LSP is therefore stable if R-parity is conserved. Cosmological arguments [^] then 
require the LSP to be neutral, and the following R-parity conserving LSP candidates are implemented 
in SUSYGEN: the lightest neutralino, the sneutrino or the gravitino. The LSP is unstable and can 
decay to SM particles if R-parity is broken, and the above cosmological arguments do not apply. In this 
case the lightest chargino, the sleptons or the third generation squarks are also good LSP candidates. 
Gluino LSPs are at present not implemented in SUSYGEN. 

Note that only a subset of the LSP candidates may be realised in constrained versions of the 
MSSM. In SUGRA models the lightest neutralino is the most likely LSP candidate, while in GMSB 
models the gravitino is the LSP. 

All sparticle decay modes relevant to LEP phenomenology are implemented with the exception 
of the three body decays of sfermions to lighter sfermions (which are only important in very specific 
regions of parameter space, see for example [33, [34|]), and the decays to gluinos. Cascade decays to 



intermediate sparticle states are also fully taken into account. 

SUSYGEN simulates the following set of R-parity conserving decay modes: 

• Three body decays of gauginos to lighter gauginos, e.g. x + ~^ ff'x°- 

• Two body decays of gauginos to sfermions, e.g. x + ~^ 

• Two body decays of gauginos to the Higgs, e.g. xt ~^ HXi- 

• Radiative two body decays of neutralinos to a photon and a lighter neutralino, e.g. xt ~~ *■ IXi- 

• Two body decays of charginos to neutralinos, e.g. x + ~^ W + x°- 

• Two body decays of sfermions to gauginos, e.g. / — » fx- 

• Two body decays to the gravitino, e.g. x° ~~ * jG- 

The R-parity violating superpotential Eq. (||) contains the 45 additional Yukawa couplings Xijk, 
Kjk> Kjk- The LLE term corresponding to X^k is antisymmetric in and we take i < j. The ODD 
term (X'-j k ) is antisymmetric in j, k, therefore j < k. To avoid fast proton decay and evade other low 
energy constraints only a subset of the 45 couplings can be non-zero. It is usually assumed that the 



A-couplings display a strong hierarchy [14], in which case it is often sufficient to consider only one 
non-zero coupling, as is done in SUSYGEN. The following set of R-parity violating decays are included 
in SUSYGEN for the three operators LLE, LQD and UDD and a single non-zero R p coupling: 

• Three body decays of gauginos to SM particles, e.g. x° ~^ for X^^ / 0. 

• Two body decays of sfermions to SM particles, e.g. l + — > l + v for X^k ^ 0. 
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3 Physics content 



This section provides an overview of the physics content of the generator. The implemented SUSY 
production mechanisms are described in Section [O], and the sparticle decays in Section [T^. The 
implementation of initial and final state radiative corrections are discussed in Section |3.3|, and sparticle 



lifetime and the hadronisation interface to JETSET are described in Sections 3.4 and 3. 



3.1 Sparticle Production 

The production mechanisms and their corresponding SUSYGEN process numbers are summarised in 
Table § in the Appendix. The diagrams for the s- and t-channel production of charginos and neutralinos 
(SUSYGEN processes 1-13) are shown in Fig. [l|. For large sneutrino and slepton masses the s-channel 
diagram dominates, while for small sneutrino (slepton) masses the t-channel contributions can be 
large, resulting in destructive (constructive) interference for chargino (neutralino) production [{nf |. 
The cross sections therefore depend on the chargino and neutralino masses and their couplings (and 
hence on Mx,M2,fi and tan/3) as well as on the selectron and electron sneutrino masses. 

Sfermions are produced in pairs (processes 14-35) through 7, Z exchange in the s-channel. The i7 e 



and e are also produced with chargino and neutralino exchange in the i-channel [17]. Fig. [2] shows the 
relevant diagrams. The sfermion cross sections depend on the sfermion masses, and on the chargino or 
the neutralino masses and couplings in the case of electron-sneutrino and selectron pair production, 
respectively. For the third generation sfermions (f, t and b) the left and right sfermion states can mix, 
and the cross sections also depend on the mixing angle 9 z |HJ . 

In GMSB models single (processes 37-39) or double gravitino (not implemented in SUSYGEN) 
production becomes accessible at LEP |20(| . The corresponding feynman diagrams are shown in Fig. |3[ 

In R-parity violating models sneutrinos can be produced singly via the LLE coupling at LEP, 
either in s-channel resonance [|35]], or in je collisions At present only the resonant production of 
sneutrinos (processes 40-45) (with non-zero couplings A121, A131) and their subsequent decays to xfP 
or Xi & are implemented. The corresponding feynman diagrams are shown in Fig. ^. 

Finally the 5 Higgs production processes (hZ, HZ, hA, HA, H + H~) are also provided through an 
interface to PYTHIA. Here the RGE improved masses and mixings of Section |Xg are used. This is 



the less developed part of the program, and we suggest that more specialised programs (e.g. see fig] ) 
are used for dedicated Higgs studies. 

In summary the 50 production processes contained in SUSYGEN are the following: 



e ' e 



XiXj 

A _ % A J 

ff 
KG 

xV 

hZ, HZ, hA, HA, H + H~ 



(21) 
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Xj (X°) 



Xt (X-) 



v e (e) 



Xj (X-) 



Figure 1: Chargino/neutralino production. 



v e (e + ) 



Xt (X 1 !) 



Figure 2: Sfermion production. 



V e (e ) 





Figure 3: Single gravitino production. 



iji 




x|(x°) 



lj (Vj) 



Xt (X-) 



e- (VJ 



Figure 4: Sneutrino resonance production. 
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3.2 Sparticle decays 



Sparticles can decay in two modes: In the R p conserving mode sparticles decay to lighter sparticles 
and SM particles (Sections 13.2.1 , 3^21 ). In the Rp violating mode sparticles can additionally decay to 



SM particles only (Sections 



3.2.1 R p Conserving Gaugino Decay Modes 

The charginos Xn an d the higher mass neutralinos xt (k=2,3,4) can decay to the three-body final 
states @: 



Xj. 



n f XmM, X° m vv 

Xk I x ± J^y, X ± 



(22) 



-J 



Fig. [5] and || show the feynman diagrams for the decays of Eq. (|2"2|). One can in general distinguish 
between two regimes: In the first the scalar masses are much larger than the gaugino mass, and the 
decays occur through the s-channel to off-shell gauge bosons and the lighter gauginos, e.g the W*x°, 
Z*x + channels in the case of the chargino decays and the W*x + , Z*X° m the case of the neutralino 
decays. In this case the different branching ratios are mostly determined from the decays of the off- 
shell W* and Z* . In the second regime the sfermions have masses close to or below the gauginos, 
and the u ; £-channels with I*, v* or q* exchange dominate and enhance the branching ratios to the 
corresponding fermions. Particular care has been taken to take into account the masses of the final 
fermions, so that e.g scenarios where one has "quasi-degenerate" chargino and neutralino masses |}7j 
give the correct branching ratios and kinematics. 
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If the sfermion mass is below the gaugino mass, two-body decays of the gauginos to the sfermions 
become dominant decay modes. This is implemented in SUSYGEN in two ways: In the first mode - 



which corresponds to the input card "TWOB FALSE" , see section |5.l| - no distinction between three 



body and two-body decays of gauginos is made. Instead the three-body calculation is used, and the 
pole at p 2 j = Mj is avoided by the inclusion of the sfermion width in the propagator^. Above threshold 
the propagator term forces the two-body kinematics, and the three-body decay is therefore equivalent 
to two two-body decays. In the second mode - corresponding to the card "TWOB TRUE" - the 
two-body gaugino decays are calculated explicitly. The following two-body decays are implemented: 




(23) 

In addition to the decay modes of Eq. (22)-(p3) the gauginos can either decay to the Higgs [pq], or 



in restricted regions of parameter space radiatively via one loop diagrams [|TJ] to x?7- The radiative 
decays are particularly important in the super symmetric limit tan(/3) —* 1; M, /i — > 0, where they are 
dominant. 




(24) 



'The propagator term D = (p'j - Mj)' 1 is replaced by D = (p'j - Ml + iGf)~ x , where Gj = FfMf. 
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In GMSB models the gravitino is the LSP, and the gauginos can decay to [21] 



X° n ^{ G 7 ,GZ,GH°,Gh,GA 

(25) 

X+ - { GW, GH + 

The coupling of the gravitino can be very small, and thus the decay Xi ~^ Gj is phenomenologically 
the most interesting mode if it occurs inside the detector. An isotropic two body decay has been 
assumed. 

3.2.2 R p Conserving Sfermion Decay Modes 

The sfermion decays to a gaugino and a fermion can also be separated into two cases. In the first case 
the lightest neutralino (x?) is the only gaugino lighter than the sfermion, and the typical decay mode 
is 

/-/Xi (26) 

with the exception of the stop, where the decay t — » t\\ is kinematically inaccessible at LEP and 
therefore proceeds through the Cabbibo suppressed mode t — ► ex?- In the second case heavier neu- 
tralinos x k = 2, 3, 4) or charginos Xn may be lighter than the sfermion, and the sfermion decays 
through a cascade of gauginos: 



(27) 




In GMSB models decays to gravitinos are also possible |21], although in practice only the decay 



of the NLSP to the gravitino is of importance. The implemented decays are: 

f Gf (28) 

The above decays cover neutralino and gravitino LSP scenarios. If the sneutrino is the LSP (or 
other sfermions are the LSP in R p violating models), sleptons and squarks can also decay via virtual 
W*,Z*,x + * an d X°* exchange to the sneutrino. At present these decay modes are not implemented 
in SUSYGEN. 

3.2.3 R p violating Gaugino Decay Modes 

SUSYGEN implements decays for a single non-zero R-parity violating coupling Ay , X[j k or \"- k , where 
k are generation indices, and the three couplings correspond to the LLE, LQD or UDD operators 
in Eq. ([]]). Neutralinos can decay to [22]: 



X° n (LLE) - { irPjl+^fu^^ilJl+^ilp^ (29) 
Xn(LQD) -> I l~Ujd k ,lfujd k ,Vidjd k ,Vidjd k (30) 
Xn(UDD) -> | Uidjd k ,Uidjd k (31) 
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and Fig. [7| show the corresponding diagrams (for example) for the LQD operator. Charginos can 
decay to |2T 



Xn(LQD) -> I ViUjd k ,lfdjd k ,lfujU k ,i>idjU k 

,UiUjd k ,UidjU k 



(32) 
(33) 
(34) 



The above decays are particularly relevant when the chargino is the LSP. Charginos will normally 
decay to the neutralino via Eq. (^) if the neutralino is the LSP, in which case the decays of Eq. (|3~^)- 
( |34| ) only dominate for large couplings A and when the exchanged sfermion mass is close to the mass 
of the chargino. 







v : 






d k " d k dj 



Figure 7: R-parity violating neutralino decays for a dominant LQD operator. 
3.2.4 R p violating Sfermion Decay Modes 

Sfermions can decay directly to SM particles through a dominant LLE coupling via 



iL 



kR 
Vi 



V 



J 



V 3 l k 
Vil k 

VilJ,Ujir 



(35) 



where In denotes a left-handed slepton of generation i. For a dominant LQD operator sfermions 
decay via 



iL 



Ujd k 
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ViL — ► djdk 

UjL -> ifdk (36) 

djL -> l^u k 

dkR -> IJUj 



dkR — > Vidj 
and similarly for a dominant U DD operator 

UiR — ► djdk 

djR -> Wj4 (37) 

Note that SUSYGEN takes into account mixing between the left-right third generation sfermion 
states. As already mentioned in Section [3.2.2 , three-body decays of sfermions to lighter sfermions are 



not included in SUSYGEN. These decays are only of importance if the produced sfermion is the NLSP 
which cannot decay directly via a specific R p violating coupling, and the LSP is another sfermion. 



3.3 Initial and Final State radiative corrections 

The initial state corrections are implemented using a factorised "radiator formula" (REMT by Kleiss) 
where exponentiation of higher orders and a Pt dependent distribution have been implemented 
They have been checked against other initial state codes and found to agree at the percent level, once 
the prescriptions of the scales e.t.c were made to agree. 

Final state radiative corrections are implemented within SUSYGEN in two ways: Firstly QED 
and QCD corrections to quark final states are treated by JETSET (through a LUSHOW call) in 
the hadronisation interface (see also section below) . QED Bremsstrahlung off leptons is implemented 
using the PHOTOS package ^(|. This allows an estimate of Bremsstrahlung effects in the leading- 
logarithmic approximation up to 0(a 2 ), including double emission of photons. 



3.4 Lifetime 

SUSYGEN calculates the lifetime of the sparticles from their decay rate, and implements secondary 
vertices for long-lived sparticles accordingly. The lifetime information is important for the R p violating 
LSP decays when the coupling strength A is small, and for the GMSB decays to the gravitino. The 
inclusion of lifetime effects can be optionally turned off. 



3.5 Hadronisation 



SUSYGEN is fully interfaced to JETSET 7.4 [24] which takes care of the hadronisation of quarks. 
Here we describe the details of the interface and the hadronisation models used by SUSYGEN. 



Two fragmentation schemes can be chosen by the card FRAG (see also section 5.1). In the 
independent fragmentation (IF) scheme the final state quarks are treated as independent particles 
with respect to each other. There is no QCD radiation between the final state quarks, and the 
hadronisation process produces relatively hard jets. Clearly this is an unphysical simplification, and 
the IF scheme should only be used to compare and evaluate the effect of gluon radiation. In the second 
scheme (the default) the final state quarks are evolved according to the Lund Parton Shower model. A 
colour string is formed between colour singlet qq states, and the two quarks are subsequently evolved 
in time/space, emitting gluon radiation between the two quarks. The jets are much softer owing to 
the gluon radiation between the quark states, which produces additional soft particles. 
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In the simplest case when there are only two final state quarks (for example in the process e + e — > 
bb — > there is no ambiguity in the colour string assignment. In the case when quarks are 

produced by colourless sparticles a colour string is formed between each of the decay products of the 
colourless states. For example in the process e + e~ —> x + X~ (x QiQ2){x°Q3Qi) a colour string is 
formed between the quarks qiq2 and another one between the quarks (73(74. 

A further complication arises when the R p violating couplings A" produce triple quark final states 
(for example in the decay x° ~^ QiQ2Q3)- The triple quark vertices are not implemented in the higher 
level user interface of JETSET, but are supported by the lower level JETSET routines. Here we use a 
prescription kindly provided by Torbjoern Sjoestrand to implement this type of vertex: The two 
quarks q\q2 are connected by a double T-junction to a diquark of flavour (qi, (72) with zero momentum 
and the third quark (73. JETSET subsequently treats the parton shower evolution of the three quarks 
correctly, but produces the error messages 

Advisory warning type 2 given after 1 LUEXEC calls: 
(LUPREP:) unphysical flavour combination 

which should be ignored. 

At present the JETSET hadronisation interface described above has two limitations: Firstly if 
squarks decay via the R p violating coupling A" to two quarks (q~i — * (72(73) the colour string connection 
between (72(73 is not supported by JETSET, and we have to resort to the independent fragmentation 
scheme instead. Secondly, because the lifetime of the decay t — > x° c is larger than the hadronisation 
scale, stops should hadronise (to stop mesons) prior to their decay. Stop meson hadronisation is not 
implemented in SUSYGEN. 
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4 Program description 



The SUSYGEN program may be used in two different ways. Firstly as a standard MC program to 
generate MC events for one particular set of input SUSY parameters. And secondly to scan SUSY 
parameter space in M2, fJ>,tan (3 and mo, where SUSYGEN calculates the masses, cross sections and 
decay branching ratios for each point, and also optionally generates weighted signal events for each 
of the specified SUSY processes. In the "SCAN" mode the calculated SUSY parameters can also be 
written to an ntuple, which can be subsequently analysed by the user. 

The user has the option to directly access the LUND common block within the program in the 
routine USER. This is particularly useful for feasibility studies in the SCAN mode, where one might 
for example estimate detection efficiencies by applying a simple set of cuts inside the above routine. 
The calculated efficiencies are stored in the ntuple, one efficiency value for each of the generated SUSY 
signals. An example USER routine which selects hadronic events with large missing Ej- is provided 
with the release version of SUSYGEN. 

In the following the structure of the program and the main subroutines are described. A de- 
scription of the common blocks is given in Appendix [D|. Inputs and Outputs of the program are 



described in sections 5^,5^. The program is divided into four stages. The first stage performs general 
initialisations: 

• reads the input cards in subroutine S CARDS using the CERN package FFREAD. The user can 
define his/her own cards in subroutine USER_SCARDS. 

• then books histograms in routine SBOOK and calls NTUPLEJNIT to define the SUSYGEN 
ntuple. 

• initialises the sparticle names in SUPART and calls the routine USERJNIT for any other user- 
dependent initialisation. 

The second stage uses the input information to calculate the masses, mixings, branching ratios 
and cross sections^. The routine SUSINI calls further subroutines: 

• SUSANA for the calculation of the masses/mixings of the sparticles. The routine SUSANA calls 
the routines: 

— SFERMION, which determines the sfermion masses either using the routine SMGUT, which 
computes the masses of the sfermions in the SUGRA model or the GMSB model, or simply by 
taking each sfermion mass from the input cards. The user can provide the mixings parameters 
of the third generation sfermions through input cards, or alternatively the sfermion mixing 
is calculated in this routine according to Eq. ( |il||l§j ) . 

— GAUGINO, which computes the gaugino masses and mixing through the diagonalisation of 
the corresponding matrices M°, M c . 

— SUBH, which computes the Higgs masses; and PINTERF, which passes the Higgs parameters 
to PYTHIA. 

— LEPLIM, which examines whether the specific point has already been excluded by LEP I 
searches. The imposition of this constraint is selected through input cards and should be 
developed by the user. 

— RPARINI, which initialises the calculations of the diagrams and decay branching ratios. 

• BRANCH for the determination of the widths and branching ratios of the sfermion and gaugino 
decays. 

5 SUSYGEN can operate in SCAN mode, where a range of input values is scanned (see input cards section). In this 
mode SUSYGEN loops over the input parameter range and recalculates the second stage parameters for each parameter 
space point. 
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— First the i? p -conservmg two body decay widths of sfermions and gauginos are calculated. 

— Routine WICONST calculates the 3-body decay constants and the function WSC calculates 
the partial decay rates. 

— Routines NTONPH and NTOXH calculate the gaugino decays to a photon or a Higgs. 

— In routine INTERF the widths of the 3 previous steps are stored in the common block 
SSMODE. The nomenclature of ISAJET 7.03 has been retained for comparison purposes. 
Then the gauge mediated and decays are calculated and stored in the same common block 
in the routines MGMDECAY and RPARDECAY, respectively. Finally the total sparticle 
widths and branching ratios are determined by summing up the partial decay rates in routine 
SSNORM. 

• The total cross sections are calculated in the routines CHARGI-HIGINIT (see table ||). The 
initial state radiation corrections are calculated in routine REMT1, and routine PROCINIT 
stores the cross sections and other information for each of the 50 processes in the common blocks 
PROC1, PROC2. 



Production 


Total cross section 


Differential cross section (da/dt) 


mechanism 


Routine name 


Routine name 


X% Xj 


CHARGI 


GENCHAR 


Y y 


PHOTI 


GENPHO 


ee 




GENSEL 


e~ie~ R 




GENSELRS 


W,TT,qq 




GENSMUS 


V e V e 




GENSNUE 


VV 




GENSNU 






GENPHOG 






RPVZINT 






RPVWINT 


hZ, HZ, H+H~, 


HIGINIT 


PYTHIA 


HA, hA 


HIGINIT + SIGHA 


PYTHIA 



Table 2: Subroutine names of the total and differential cross-sections. 



In the third stage "weighted" MC events are generated, i.e. one of the selected production processes 
is chosen with a probability proportional to its relative cross section. The routine SUSEVE then loads 
the corresponding variables from the commons PROCl,PROC2 and performs the following tasks: 

• Picks a value for cos 9 distributed according to the differential cross section {da/dt) for the process 
(see also table ||) and creates the 4-vectors for the final states of the hard process. 

• The decay channel is determined in the routine DECABR from the branching ratios of the spar- 
ticles, and the two-body or three-body decay 4-vectors are calculated in the routines SMBOD2 
and SMBOD3, respectively. In SMBOD3 the differential decay rate dT/dsdt is sampled in s and 
t, and the angular distributions of the decay products determined accordingly. 

• The decay products from the previous step are presented again to DECABR and decayed, until 
no more unstable SUSY particles are left. 

• The routines REMT3 and FINRAD_PHOTOS perform the adjustments due to initial and final 
state radiation. 
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• The above 4-vectors are interfaced to LUND in the subroutine SFRAGMENT where they frag- 
ment, hadronise and decay. 

While still in the third stage, inside the event loop, the routines USER_WRITE and USER are 
called where the user has the options to write an event out and/or analyse it directly by accessing the 
LUND common LUJETS. 

At the end of the event loop - while still inside the SCAN loop - the calculated parameters (masses, 
BRs, cross-sections) and the rough detection efficiencies calculated in the routine USER are stored in 
an ntuple by the routine NTUPLE_FILL. 

In the fourth stage the routines SUSEND and USER_END perform general and user-specific closing 
functions. 
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5 Setting up and running SUSYGEN 



The program is supplied as a single Fortran file which can be compiled using standard FORTRAN77. 
The release version may be obtained from the CPC Program Library, and the most recent version can 
be obtained from http:/ /lyohp5. in2p3.fr/delphi/katsan/susygen. html. The program has to be linked 



to the standard CERNLIBs, the JETSET7.4 library, PYTHIA and PHOTOS. This may be done by 
the following set of commands on UNIX: 

f77 -g -w -static -c susygen.f 

f77 -o susygen.exe susygen.o /cern/pro/lib/libphotos . a \ 
/cern/pro/lib/libjetset74 . a 'cernlib mathlib packlib shift' 

and on VMS: 

$ f or/optimize/nodebug susygen 

$ cernlib j etset74, photos , packlib, genlib,kernlib 
$ link/nodeb susygen, ' lib$ ' 

To run the program one has to give appropriate input cards, which in the example below are 
assumed to be in the file susygen. cards. To run SUSYGEN on UNIX: 

susygen.exe < susygen. cards > susygen.log 

and on VMS: 

$ define/user sys$input susygen. cards 
$ define/user sys$output susygen.log 
$ run/nodeb susygen 

The Inputs and Outputs are discussed in the following. 



5.1 Input 

The input to SUSYGEN consists of a single text file, the cards file, described below. Alternatively 
a very advanced graphics X-interface, Xsusy, is also available for UNIX machines. The Xsusy inter- 
face allows input from pop-up menus rather than a cards file, and is extremely user friendly. The 
SUSYGEN output may be analysed with PAW and/or emacs within Xsusy. Fig. || shows a typical 
Xsusy-session. Xsusy is not part of the standard CPC distribution of SUSYGEN, but may be ob- 
tained from http: / /lyohp5.in2p3.fr/delphi/katsan/susygen.html . 

In the absence of any of the cards the default values given below are chosen. 



5.1.1 Set masses and mixings 

MODES 1 MODES determines which input parameters will be used to determine the supersym- 
metric masses and mixings: 

• MODES=l the gaugino/higgsino and sfermion spectrum is calculated using the values M, fj,, 
mo, tan/?, and A as input for the case of SUGRA models, and M, fi, tan/? and A for GMSB 
models. 

• MODES=2 the gaugino/higgsino spectrum is calculated using the values M, /u, mo and tan/? 
as input, while the the sfermion masses are taken from the input cards below. 
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• M0DES=3 the gaugino/higgsino spectrum is calculated using the values M, fi, mo and tan (3, 
but the "gaugino mass relationship" does not relate M\ to M2 anymore, their relationship is 
set by the card RS below. The sfermion values are given according to the hierarchy proposed 
for MODES=l. 

• MODES=4 the gaugino/higgsino spectrum is calculated using the values M, fi, mo and tan (3, 
but the "gaugino mass relationship" does not relate Mi to M2 anymore, their relationship is 
set by the card RS below, while the sfermion values are taken from the input cards below. 

• MODES=5 TO BE USED WITH CARE. M below is the requested mass of the lightest 
neutralino. A corresponding M2 may not exist for the given set of parameters, or very often 
double solutions are possible (e.g for positive fi and low tan (3). 

M 90. (M 2 the SU(2) gaugino mass in GeV,at the EW scale) 

mu 90. (/x in GeV, at the EW scale) 

mO 90. (mo in GeV, a common scalar mass at the GUT scale) 

tanb 4.0 (tan/3 = vijv\, at the weak scale) 

A 0. 0. 0. (At, A;,, A T soft SUSY breaking trilinear couplings) 

RS 1. RS is the gauge mass unification breaking scale defined as M\ = M ^^^^ RS) 
MA 500. (The mass of the A Higgs in GeV) 
MIX If set to 1, permits sfermion mixing. 

C ARE RELEVANT FOR MODES 2 and 4 

MSQUARK 1000. (the mass of q, m~) 
MLSTOP 1000. (the mass of t L ,m f ) 
MRSTOP 1000. (the mass of t R , mr ) 
MLSEL 1000. (the mass of e~ L ,m e ~ L ) 
MRSEL 1000. (the mass of e~ R ,m e ~ R ) 
MSNU 1000. (the mass of 0, m~) 

MASMX 0. 0. 0. 0. 0. 0. When different from 0. sets the mass of t\, bi, f\ £2, 62, "?2 respectively. 
PHIMX 0. 0. 0. sets mixing angle between t\-t2, b\-bi and f\-f<i- 

5.1.2 Set production processes 

ZINO FALSE (all pairs of x°) 
WINO FALSE (all pairs of x ± ) 
SELECTRON FALSE (all pairs of e) 
SMUON FALSE (all pairs of fi) 
STAU FALSE (all pairs of f) 
SNU FALSE (all pairs of v) 
SSTOP FALSE (all pairs of t) 
SBOTTOM FALSE (all pairs of 6) 
SQUARK FALSE (all pairs of q other than t, b) 
HIGGS FALSE (all higgses a la PYTHIA) 
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GRAVITINO FALSE (Gx°, when MGM is true) 

RSINGLE FALSE (ux° and or tx ± through v resonance when RPARITY is FALSE) 

PROCSEL 0000000000 

Possibility to choose to produce any 10 processes. The first word is the number of processes 
to be selected and the next ten give the process number (see process table in the text). The 
corresponding generic card is turned TRUE. That is e.g ZINO is set to TRUE when one sets 
PROCSEL 1 2, asking to produce only the process x\x%- 

DECSEL 11111111111111111111111 

Determines the modes of decay of virtual Z and W. The first 4 determine whether the virtual Z* 
can decay to uu, dd, uu, e + e~ and the next two the decays of W* to ud or e + u. The next 12 
determine the decays to the two other families. The last 5 determine whether one wants to turn 
on/off the decays to xa, Xih,XiH, Xi A , xfH T . 

5.1.3 Gauge mediated SUSY models 

MGM FALSE Switches on Production (%£?) of gravitinos and the decay modes to the gravitino 
in models of gauge mediated supersymmetry breaking. 

LSUSY 500. Supersymmetry Breaking Scale y/F (in GeV) (gauge mediated SUSY models) 

5.1.4 R p violating SUSY models 

RPARITY TRUE Is R p conserved? 

INDIC o i j k R p violating coupling indices of the single dominant coupling \jk- The operators 
o=l,2,3 correspond to the Rp violating operators LLE, LQD, UUD respectively. 

LAMDA 0.1 Strength of the R p violating coupling Xyk- 

TWOB FALSE Switch on/off explicit calculation of the gaugino two-body decays. The TWOB 
is automatically set to true if R p is violated. 

5.1.5 Generation conditions 

ECM 183. (Center of mass energy). 

SUSEVENTS 100. (Total number of events to be produced. If more than one processes are 
possible, they will be produced in percentages corresponding to their respective cross-sections). 

ISR 1 (Permit or not (0) initial state radiation). 

FSR 1 (Permit or not (0) final state radiation). 

GENER (Generate events(l) or not (0)). 

LIFE TRUE (Include effects of finite lifetime of sparticles). 

FRAG 1 (Fragmentation scheme. Normal string type fragmentation (1) or independent fragmen- 
tation (2)). 

LCUT 200. (Decay length cut in cm, thereafter particles considered stable). 



22 



5.1.6 Control switches 

DEBUG (Print additional debug information. The higher the value, the more information is 
printed (0 <Debug< 31)). 

LUWRIT FALSE (Write out events to FOR012.dat or not) . 

SCAN FALSE (Scan the values in M and fj,). 

ISCAN (Number of scan values in M2, fJ>, mo and tan/3). There is also the following 
option: setting the relevant ISCAN value to -1, the program takes the first corresponding value 
VSCx as a starting value, the second as ending value and the third as the step. 



VSC1 





(List 


of scan values in M2, up to 10, or start, end, step values) 


VSC2 





(List 


of scan values in /i, up to 10, or start, end, step values) 


VSC3 





(List 


of scan values in mo, up to 10, or start, end, step values) 


VSC4 





(List 


of scan values in tan f3, up to 10, or start, end, step values) 


LEPI 


FALSE 


(Impose LEPI constraints or not). 



5.2 Output 

Four output files are produced by SUSYGEN: 

• The file susygen.dat, which lists the sparticle masses, cross sections, branching ratios and decay 
modes. 

• The file susygen.his. The ntuple with ID 200 in this file contains the calculated masses, cross 
sections and BRs. 

• The fortran file UNIT 12, e.g. "fort. 12" on UNIX or "FOR012.DAT" on VMS, which contains 
the LUND common of the generated events. See routine SXWRLU in the program for more 
information on the output format of this file. 

• The file susygen.log, which contains warning and error messages and a LUND listing of the 
generated events (if the debug level is set to "DEBUG n", where n > 0). 

Ntuple 200 is a column-wise Ntuple, and contains one entry per SCAN loop (SUSYGEN can scan 
parameter space in /i, M2, tan/3, m_o). The ntuple variables are summarised in Tables |Tll-p2| in the 
Appendix. 



23 



Xsusy 



Steering Panel 



J 



Input Parameters 



u 



Production 



JJ. 



X Susy gen 2. 2 
Input parameters 



Compile Susygen 



Show Results 



Susygen Manuel 



Heip 



Quit 



Steering Panel 
Default Values 



old.csh 



Set SUSY model 
♦ SUGRA or v GMSB 

Set SUSY Spectrum 



R-pariiy 



Set production processes 



Set Run conditions 



Expert 



CreateScript hew.csh 



Help | Quit | 
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i; zone 2 2 

': title 'SUSYGEN 2.2 / SIJGRfi m0=80, tan[b]=1.41, I 
i: 2dhisto 10 'mchi' 41 -200 200 20 10 200 
l: nt/loop 200 f.ftfx01,f»2,fnu,10.) 

h/plo 10 surf 
.; atitle '[»] (GeV)' 'M2 (GeV)' 'M[h] (GeV)' 
:: nt/plo 200,(100.*(brxplb+brxplc+brxpld)/(brxpla+t 
:)»fn2 friu=-200 ! ! ! L 
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Figure 8: The graphical user interface Xsusy. 
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6 Conclusions and Future plans 



SUSYGEN 2.2 is a fast and versatile Monte-Carlo generator that has profited from the collective 
experience of the four LEP collaborations. Its current precision is sufficient for detailed studies at 
LEP. It permits the study of gravity and gauge mediated models of SUSY breaking in models in which 
R-parity is either conserved or violated. A total of 50 production processes are included, which are the 
most important processes in the models studied above. SUSYGEN is extremely user friendly and the 
code is transparent. It is maintained in CMZ, and private developments can easily be incorporated in 
the main body of the generator. 

There are a few things the authors would like to incorporate in the next version: 

• Extend SUSYGEN to hadron and ep colliders. 

• Treatment of helicity correlations. This development is well under way and should be imple- 
mented in a few months. 

• Include a light gluino scenario. 

• Include single squark/sneutrino production in fi p processes where a radiated photon from one 
of the incoming leptons interacts with the opposite lepton through a "resolved" quark or lepton 
component. 

• Include three body decays of sfermions to lighter sfermions. 

• First order corrections to SUSY masses and cross-sections. This development could be important 



in particularly difficult areas of the parameter space (e.g. tan/3 = 1 see [45]), where corrections 
of the oder 20% are possible. 

Other more long-range plans include: 

• A more precise Higgs treatment, which also includes Higgs decays to SUSY particles. 

• Incorporation of cosmological constraints and other constraining measurements (e.g. b —* S7 
constraints) at an informative level. 
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APPENDIX 



A Codes for sparticles and processes 

Table || lists the LUND KF codes which SUSYGEN assigns the sparticles in the LUND record. Note 
that SusyProd stands for the production process of the sparticles, upper case G is the Gravitino and 
lower case g is the gluino. Table ||| lists the SUSY production process code conventions. 



LUND KF id 


Sparticle 


LUND KF id 


Sparticle 


41 


d L 


61 


b R 


42 


UL 


62 


in 


43 


SL 


63 


h 


44 


CL 


64 


h 


45 


h 


65 


h 


46 


h 


66 




47 


da 


67 


h 


48 


UR 


68 




49 


SR 


69 


G 


50 


CR 


70 


9 


51 


el 


71 


X? 


52 


V e 


72 




53 


P>L 


73 


x° 3 


54 


v» 


74 


xl 


55 


f L 


75 


xt 


56 


V T 


76 


xi 


57 


Z~R 


77 


Xi 


58 


&R 


78 


Xi 


59 


f R 


79 


SusyProd 



Table 3: LUND KF labels and particle codes (K(I,2) in the lund common). 
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Process id 


Process 


Process id 


Process 


1 


x u ixi 


26 


Mi 


2 


xU 


27 


hh 


3 


X°2X°2 


28 


ulu l 


4 


xUi 


29 


URUR 


5 


xixi 


30 


dhdh 


6 


xlxl 


31 


ArAr 


7 


x?x2 


32 


CLCL 


8 


xlxl 


33 


crcr 


9 


xlxl 


34 


s~ls~l 


10 


xixi 


35 


S~RS~R 


11 


xixi 


36 


X°rG 


12 


X1X2 


37 


X°2G 


13 


X2X2 


38 


XIG 


14 




39 


X\G 


15 




40 


Xx u 


16 


i7 T u T 


41 


X% v 


17 


e~L(iL 


42 


X% v 


18 


e~Re~R 


43 


x\v 


19 


e~Le R 


44 


x\\- 


20 




45 




21 


mm 


46 


hZ 


22 


TLT L 


47 


HZ 


23 


TRTR 


48 


hA 


24 


hh 


49 


HA 


25 


b 2 b 2 


50 


H + H~ 



Table 4: SUSY production process codes. 

B Gaugino Diagonalisation 

B.l Mass matrices 

The neutralino mass matrix in the base 

$ = (-#y, -iZ, cos (3Hf - sin 0H$, sin 0H$ + cos 0H^ (38) 

is given in section 

|3 The mass eigenstates can be computed by diagonalising Eq. (^). The diago- 
nalisation is done numerically in SUSYGEN. In some calculations (R p violating processes and gauge 
mediated processes) we use the base ]42[ |: 

i>k = {-^B'-^Wa^Hl^m) (39) 
In this base the diagonalised neutralino coupling matrix N-j is defined in terms of the matrix Nij as 

N[j = N 1:j cos(6 w ) - N 2j sm(6 w ) 
N' 2j = N lj svn(e w ) + N 2j cos{e w ) 
N'v = N 3j cos(/3) + N 4j sin(/5) 

N' 4j = -N 3j sm(f3) + N 4j cos(l3) (40) 
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The chargino mass matrix is also given in section 2A. Here the chargino masses can be computed 
analytically: 

Ki,xt = \[M! + fi 2 + 2m 2 w (41) 



TV(Af| - /^ 2 ) 2 + 4m^cos2 2 /3 + 4m^(M| + fi 2 + 2M 2 /xsin2/3)] 
The chargino mixing matrices are given by 



Ua = V21 = ^J 1 + Mi-^-^^ {i2) 

Un = _ gu= »^MEZ ^~^ ,43) 

K. = ^.A^Z^w (44) 

f 22 = v» = ^^/^EZ+pi-iZ (45) 



where W is given by 



W = J (Ml + fx 2 + 2m 2 w ) 2 - 4(M 2 /i - sin2/3) 2 (46) 



and the sign factors 9i are shown in table ||. The chargino and neutralino mass eigenvalues , m x o 
can have positive or negative signs % = ±l,r] n = ±1, while the physical masses are defined positively, 
e-g- M + = m +r] k . 

X-k X-k 





tan (3 > 1 


tan (3 < 1 


01 


1 


£b 


e 2 


£B 


1 




£A 


1 


04 


1 





Table 5: Sign factors 6i, where e yi =sign(M2 sin /? + /icos/3) and £e=sign(M2 cos (3 + //sin/3). 



31 



C Formulae 



In the following we list the formulae which are used in SUSYGEN for the cross sections and the decay 
rates. The listing is given for completeness and easier understanding of the code, together with the 
original references. 

C.l Cross sections 

C.l.l Neutralino pair production (processes 1-10) 

Differential neutralino cross section (function GENPHO) fl7|| ; 

da da da. da. , . 

~dt = H lz+ ^ + H lzi {1) 



da g 



D z {s)\ I \O il L \\\L e \ i + \R e \ 2 ) (48) 



[{Mf - t){Mj -t) + (Mf - u){Mf - u) - 2ri i r lj M i M j s] 



da. g 



{|//f| 2 |/^ 2 [l^(i)| 2 (^ - t)(Mf - t) (49) 



+ \D e ~ L (u)\ 2 {M* - u)(Mf -u)- 2Re(D e ~ L (t)D* e ~ L (u))7 liVj M i M jS ] 
Mfi?\V$\ 2 l\D e ~ R (t)\\M? -t)(Mf -t) 

+ \D e ~ R (u)\\Mi - u)(M? - u) - 2Re(D e ~ R (t)D* e ~ R ( U )) VlVj M t M,s]} 

{LJ^f^D^mM? - t)(Mf - t) - mVjMiMjs) 
+D e ~ L (u)({Mi - u)(Mf -u)- mVj MiM jS )] 
-R e fXf§[D e ~ R (t)((Mf - t)(Mf -t)- m^MiMjs) 
+D e ~ R {u){(Mf - u)(Mf -u)- mrijMiMjs)]} 

where the Standard model couplings are; 



9 
Lf 



e 

sin #w 

T 3 - Q sin 2 W 
— Qsin 2 #w 



(51) 
(52) 

(53) 
(54) 



Mi and r\i are neutralino masses and diagonalisation phases, the couplings of the neutralino to 
sfermions are; 
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fk = -V2[^-N t2 --^ Nll ] 
cos t/\v sm "w 



the neutralino couplings to Z are; 



<V = - 2 ^ N h " ^4iV* 4 ) cos 2/3 - - (N l3 N* A + ^ 4 7V* 3 ) sin 2/? 
"R = _ "L* 

and the propagators are; 



Dz(s) = (s - m| + imwTz) 1 
D~ p (x) = (x-m 2 ~ )~ 1 



The integrated cross section (function PHOTI) is; 



2 8ij . . 
o-tot = — ~ — (o-z + o-g + a ze) 



4 2 

CTZ = 47rco 7 s 4 g w |jPz(s)|2 7i |Q ^ |2(|Le|2 + l^l 2 )^ + Y ~ ^ M * M il 



4 

9 1 



9 



+l/5l 2 l/5l 2 [^^^ + 2 + ^£(i-2d B - 5« 

7 sdi. - q 2 2q sd R d R *= 

V s 



— = i? e ( J D z ( S ))o;; L 



87TCOS Oyf^/s 

Wi/fl ^ - ^ - f - ™> M ' M ' H-^\ + 2,1 - dL)] 

-us® ** 1 - ' M1 : ~ "^ h^i + 2d - *,)]} 



the phase space factors are; 
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d L ,R = ±- s (s + 2m 2 ~ LR -M?-Mf) 



Ei = ^q 2 +Mf 
where q is the CM momentum of the neutralinos. 

C.1.2 Chargino pair production (processes 11-13) 

The differential chargino cross section (function GENCHAR) is |l7|; 



da da da da da da da 

It = It ' 7 + H |z dt lp + dt ,7Z dt w It lzp 

^It = ^[(Mi-u)(Mf-u) + (Mi-t)(Mf-t) + 2M i M j s] 



+(M? - t){M] - t)] + A{L 2 e + F^OifOifrHVjMiMjB 

-(Ll - Rl){\0'^ - \0[f\ 2 )\{M 2 - u){M] - u) - {Mf - t)(Mf - *)]} 

-\V iX \*\V jx \ 2 {M? -t){M? -t) 



da g*\D„(t) \ 2 Mr ,2, 



+(M 2 - t)(M 2 -t) + 2MiM jS ] - (L e - R e )(d± - 0[f) 
[(M 2 - u)(M 2 -u)- (M 2 - t)(M 2 - *)]} 



where the chargino couplings to the Z are; 



0$ = -VaV^ -IvaVfr + Sy sin 2 <7 AV 
1 



0<f = -^^--C/^ + ^sin^w 
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and; 



D D {t) = (t-ml)- 1 

The total cross section is; 



J^Dz(s)|^ , L 2 2 2 2 

(£^- + U + 2(L2 + Rl)0' i j'0' i fr]i rjj Mi Mj ] 



g 4 l^il 2 l^il 2 g f ^- + g 2 -g^t | V 

32vr^mi [ a| - b\ + b\ 



i 



2bi b L a L - b L 



(£7^ + £ + MiMj) 



e 2 g 2 \Vn\ 2 S ij | ql + &l | 

<V = — + — 



® zv = - 167TS 9 cos 2 ^ iVjiRe{Dz(8))L e (0'^h + C^f r^M^lnl^. 



where the phase space factors are; 



2m£ + s-M? -Mf 

°i - ^2 



2ro£ 



m~ 
v 



h = 2q^/s- 2q — + {EiEj + q -f - q^/s— )ln\ — 

bL bi b L a L - b L 



35 



C.1.3 Sneutrino production (processes 14-16) 

The differential cross section (function GENSNUE) for v e is [|17|]; 



fe=l 

2L 2 



ite(£z(s))X!l^d 2£> k(*)] 



cos 2 6> w ,. ] 

where the t channel chargino propagators are; 



D k (t) = (t-Ml)- 1 
For v other than v e (function GENSNU) only the first term is retained. 



C.1.4 Selectron production (processes 17-19) 

The differential cross section (function GENSEL) is fl7|| ; 



da da da da da da da 

H = dt ' 7 + ~dt ]z It '* + It ,7Z ~dt 'tx + dt lz x 



We obtain ; 

(i) e + + e~ — > + e{ 



da e 4 4 



^|z = i 67rg2 f os 4 % ( ut ~ rn^lDzWWLe cos 2 + R e sin 2 </>) 2 (£ 2 + R 2 e ) 



= ^{(^-<)[^^(E D *(<)l/4l a ) 2 + Bin 4 0(E D *(<)l/5l 2 ) a ] 

4 

+2 S sin 2 </>cos 2 ^ M k M mm D k (OA (t)f£fekfSfel} 
k,l=l 

-j^Uz = ■^ r ^^^(ut - m^)(L e cos 2 (f> + R e sin 2 4>)(L e + R e )Re(D z (s)) 



da p 2 n 2 4 
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da 4 

iM = 1^ w e w {ut ~ m ^ )(Le cos2 * + ^ sin2 ^ e(IMs)) 

£ A(t)(L e |/i| 2 cos 2 + iiel/iH 2 sin 2 0) 
fe=l 

where the i channel neutralino propagators are; 

D k (t) = {t-Mir 1 

and (/) is the L/R mixing angle; 

ei = cos + e# sin 0, 
&2 = 6r cos — ei, sin 

(ii) e+ + e~ — > + (function GENSELR) 

dcr dcr da 

~dt ,7= rft" ,7Z= dt w = 



= 167rs2 f QS 4 9w ("* ~ «) sin 2 ^cos 2 4>\D z (s)\\L e - R e f{L\ + i? 2 ) 



= ^{K--l i ml 2 )c O s 2 0sin 2 0[(E^(t)|/i| 2 ) 2 + 

fc=i 

(^L> fe (t)|/«| 2 ) 2 ]+ S (sin 4 ^ + cos 4 0) ^ M k M m D k {t)Dmf? k fekfeifei} 
fc=l fc,z=l 



dcr £i 4 



= 167r/cOS 2 % ^ ~ ^i^^ CQS2 ^ Sin2 ^ - R e)R<D Z ( S )) 
4 



E^*(*)(^e|/4| 2 --Re|/5| 2 ) 
fc=l 

C.1.5 Other charged sfermion cross sections (processes 20-35) 

The differential cross section (function GENSMU) is; 



da da . da . da . 

~dt = ~dt ' 7 + ~dt lz + ~dt lyZ 
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(103) 



d 4 

-£\z = W7T J cos 4 0w ~ m 4 j)\Dz(s)\\L f cos 2 cp + R f sin 2 4>f{L 2 e + R 2 ) 

(104) 



Q 2 e 2 (5 3 
8tts 6 



Where /3 = H 

For q the above cross sections are multiplied by a QCD factor R; 



C.1.6 Neutralino-Gravitino production (processes 36-39) 



The cross section is given by [20 



da 1 F(s,t,u) 



dt 16vrs 2 A| 



(105) 



da Q 2 g 2 

— \jz = - — n T7^~( ut - m % )( L f cos2 4> + Rf sin 2 (f>)(L e + R e )Re(D z (s)) 

dt 87rs d cos &w Ji 

The total cross section (function GENSMUS) is; 

a = <r 7 + cr z + <r 7 ^ (106) 



(107) 



^ = i« gS 4 fl |^(^)| 2 (%co S 2 + % S m 2 0) 2 (^ + j; 2 )^ (108) 

167TCOS % 



Q2 2 o3 

a 7Z = (L f cos 2 <P + R f sin 2 <P)(L e + R e )Re(D z (s))^- (109) 

8-7TCOS #\y O 



4ry 7r 2 7r 2 
A = 3(1 + -^(^ - (1 + 0)( T - 3))) (110) 



(111) 



A| = 6(Mplanck'>tt-g) 2 is the supersymmetry breaking scale, as a function of the Planck mass 
and the mass of the gravitino. 

F(s, t, u) = F 77 + F tt + F uu + F ZZ + F~ ft + F 7 „ + F zt + F Zu + F lZ , (112) 
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where 



and: 



„ 2s(s-ml)(t 2 +u 2 ) , 
F 77 = (N iX e) 2 ^ (113) 



s 

t 2 (m 2 x -t)(-t) t 2 {m 2 x -t){- t) 

u 2 {m 2 — u)(—u) „ u 2 (m 2 — u)(—u) 



( X R) 2 n X 2 2 + ( X lV h * 2 2 ( 114 ) 



Fun = {XrY x 2 2 ' +{X L ) 2 * 2 2 (115) 
{u-mj R ) z (u — m~ L y 

Fzz - r 2 ^;J —2— ( a - + (rzM z) 2 (116) 

^ = (^eX,)i^-(^ L )i-^ (117) 
s(t-mf R ) s(t-mf L ) 



F 1U = (N a eX R ) \ - (N a eX L ) \^, ( (118) 



Fzt = -[N i2 ReX R 



(-u)(2su 2 ) , (-u)(2su 2 ) 

— 2— - (iViieAzJ— 2— 

s(u — m~ R ) s(u — m~ L ) 

g \ {-t){2st 2 )(s - M 2 ) 



cos6 w J [{s-M 2 ) 2 + {T z M z ) 2 ]{t-ml R , 

| V 7 v 9 \ (-t)(2st 2 )(s-M 2 ) 
T ( Nl2 LeXL ^) [( S -M 2 ) 2 + (T z M z ) 2 ](t-ml) (119) 

r/ii __ Xt2R<XR 9 ^ (— u)(2su 2 )(s — M§) 



cos^y [(s - Mf) 2 + (T z M z ) 2 ](u - m\ R ) 



cos 9 W J [(s -M 2 ) 2 + (T Z M Z ) 2 } (u - m 2 - ) 



+ lN i2 L e X L —S—) Tl ,^ (2 ^'|!-5 ) (120) 



(121) 



L ' 

V - -of AT N 9 ^ ^ ± 2S(S ~ m ' )(t2 + U2){S ~ M Z 

JZ ~ Z V^ ae cos9 w ) 2 s[(s-M 2 ) 2 + (T z M z ) 2 } 



X R = N lie -N i2 -^- (122) 
cos Ow 

X L = -Na e - N i2 — ^— L e (123) 

cos Ow 

C.1.7 S-channel Sneutrino resonance cross sections (processes 40-45) 

The R p violating single sparticle production processes 40-45 are characterised by a sneutrino s-channel 
resonance for a non-zero coupling A121 or A131. The s-channel sneutrino process e + e~ — * Uj has three 
principal decay modes: 

e + e~ -> Uj -» e + e~ (124) 
e+e" -» ^ -> X % (125) 
e + e - ^ ^ ^ (126) 

but only the single chargino/neutralino processes (|125j),(|i~2"6|) are implemented in SUSYGEN. 
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C.1.8 Single Chargino Production 

Ignoring contributions to the vertices of the MSSM from mass terms, we have two channels present 
(s and t) for chargino production [ 23 1 : 

1 ( V m i\ijie \ 2 \ C s C t _ 

°e+ e --^x± 64tts 2 V sin% J 1 (s - m 2 -.) 2 + rgm 2 -. (t - m? ) 2 



2 



C st (s-m^) 



[( s _ m 2 )2 +r 2 m 2 ](t _ m 2 J 



(127) 



where 



C s = s(s-m~±), C t = t(t-m-±) 

C st = s{s-m\ ± ) + t{t-m\ ± ) -u{u-m\ ± ) (128) 
Here s, t and it are the Mandelstam variables 

t = {e~-x ± )\ u = (e + -x ± )\ (129) 
and the width of the sneutrino is given by 

= ri + r 2 (130) 

A 2 

where the R p violating sneutrino decay rate T\ = uj — > e l e k = -j^ m Uj] and the R p conserving decay 
rate T 2 = vj -> Xm l J ■ 

C.1.9 Single Neutralino Production 

Ignoring mass terms the cross-section for this process is given by p3|: 



1 ,2 ta2 C s , T r 2 Ct T -o C« 



a e+e--,yx°m ~ i^2 A lil 1 ^ f e _™2 A2 T r 2^2 + ^2 (+ _ m 2 \2 + ^ 



32tts 2 131 \ 1 (s - ml-.) 2 + r?m| 2 (t-m?J 2 d (« - m?J 2 
C sf (a-m|.) ^ C sw (s-m|.) 

_ ' 2 [(s-m^)2 +r 2 m y (t _ m |j " 1 3 [{s _ m 2^ +T 2 m 2- [{u _ m 2 R) 

a 



tu 



"W ( i_ m 2j (n _ m 2j 



(131) 



where 



I cos t>w 

cos t>w V ^ 
3 3m+ cos6 w m2 



(132) 
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Here the functions C s , Ct, and C s t are given as before, just replacing the chargino mass m^± by 
the neutralino mass m^o. The new functions C u , C su and C u t are given by 

C u = u(u-rrv~ ) 

C su = s(s — m| ) - t(t — m| ) +u(u — m| ) 

C tu = -s(s — m| ) + t(t-m|o) + it(n-m|o) (133) 
And the i? p conserving decay rate T 2 = £j — ► Xm^j- 

C.2 Decay Rates 

C.3 Gaugino-Higgsino three body decay rates 

The general 3-body decay from a chargino/neutralino to a chargino/neutralino and 2 fermions can be 
parametrised as [o]; 

r(x<— >X fc + //) = QO ."A Ar3 / + W t + W„ + WW + W st + W su ) 

327TS1H Vw-M? J 

(134) 



Ws = I- ^ n 2 {MMf - m~ Mi) (135) 

(s - M^) 2 + G s z 

+5 S (M 2 - n)(n - M 2 ) + 2C, Wfc M i M fc f} 
w L (Mf-t)(t-M|) a (jg-jKg-jg) 

L (M 2 -n)(n-M 2 ) fl (M 2 -^-M 2 ) 
" " U {u-MlY + G^ +K (u-MlY + G^ (137) 

W tu = 2^^M fcS - [(f _ ^ ^ (138) 

+2^» n MM* (i - M 2 R )(u - M 2 R ) + G R 2 
+ 2A M M t M k s [{ ._ m2r)2 - Gr2][( _ _ M , )2 - Gr2] 

^ = 2^(M 2 -t)(t-M 2 )- (t ""f l)(S ";5 ) + ^f S 21 (139) 

st [(t-il^) 2 + G L 2 ][(«-M?) 2 + G B 2 ] 

+2vl R (M 2 — t) (i - Mi) ~ M 2 R )(s - M?) + G R G S 

+ A 1 Kt Mfe) [(t-M 2 )2 + G R 2 ][( S --M|)2 + Gs 2 ] 

, „„» „ M M , (t - M 2 R )(s - M 2 ) + G R G S 
+ 2B M M iMk s - _ m2r)2 + Gr2] [( _ _ m|)2 + G , 2] 

W - 2A L (M 2 u)(u M 2 ) (u-M 2 L )(s-M 2 ) + G L G s 

.oj,L nn MM, (u - M 2 )(s - M 2 ) + G h G s 
+ 2B M M lMk s - _ m2)2 - Gl2][( _ _ _ - Gs2] 
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+2A? u (M?-u)(u-M 2 k ) 



-2B*rHT} k MiM k 8 



(u - M 2 R )(s - M 2 ) + G R G S 
{{u - MlY + G R 2 ][( S " - MIY + G 2 ] 
(u - M 2 R )(s - M 2 ) + G R G S 
[{u-M 2 R ) 2 + G K 2 ][{s-M 2 ) 2 + G 2 } 



where the factors G { are; G s = T B M S , G L = T L M L , G K = T R M R , G v = T U M 9 , G k 
G\ = T\M\ and the chargino-neutralino couplings are; 



r k M k , 



L 

o R 



v* 

- (N u cos - N i3 sin 0)^% + (Nn sin 6 W + N i3 cos %)V£ 
(iV* 4 sin + N* 3 cos P)^ + (iV*! sin 9 W + N* 3 cos 9 W ) Uj i 



(141) 
(142) 



and the tables [| |8| and |9| contain the definitions of the A, B and C parameters. 



Xi 



l+l- 



x k vv 
cos 4 6*w 
cos 4 6>w 
cos 4 6*w 



(f/fci^i) 2 






cos 2 9~w 



COS' 



A, 



A L 
A L 

A L 



B L 



B 



2(0£) 2 

(/4^i) 2 

-y/2tf k U ix Ok 



cos 4 #w 

cos 4 6*w 

cos 4 0\y 
(Hi^i) 2 





cos 2 6*w 
cos 2 #w 





Table 6: Coefficients for leptonic chargino decays 



C.4 Gaugino two body and radiative decay rates 

The above three-body branching ratio formulas include the width of the exchanged particles in the 
propagator terms, and therefore smoothly converge to a two-body formula when above threshold. For 
example if the sneutrino is lighter than the chargino and the neutralino is the LSP, the three-body 
decay rate T{\ + ~^ l +l/ X°) 1S the same as the two-body decay rate T(x + ~^ Furthermore the 

kinematics of the three-body decay is equivalent to the two-body decay kinematics (once the sneutrino 
decay v — > ux° has been included). The procedure has been numerically checked with the explicit 
2-body calculations, and the branching ratios, partial decay widths and kinematics agree well. 

Sometimes it is useful to calculate the gaugino two-body decays explicitly (for example for R p 
violating models). The gaugino decay rate to a lighter sfermion is given by: 

r(x -►//') = c f l6 ^ m3 \J {m 2 x - (mj + m f ) 2 )(m 2 x - (mj - m f ) 2 ) 
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xl ud 



^k 

cos 4 #w 
l2 ^o'K?+{RuO' k \? 

cos 4 #w 
12 (Ll+Rl)0't0' k f 

cos 4 6*w 



3(f/ fc i^i) 2 





Q U k iUgO' k fL u 
cos 2 #w 

77i 

cos flw 



A s 
B s 

C s 

A\ 
A L 



Bit 
A L 



B 



6(0^) 2 
3(/4^i) 2 

3(/i^i) 2 

ZfukfkVaU* 

3V2fgUaO* 
-3y/2fgU a Og 



12 (L d O k ^+(R d O k f)^ 

cos 4 #w 

1 o (^oL?) 2 +(fl d O^) 2 
iz 

cos % 

12 (^+*2x&q£ 

cos 4 #w 
3(F fe i^i) 2 



Q VkiViiO'£L d 

cos 2 0\y 
cos 2 #w 






Table 7: Coefficients for hadronic chargino decays 



Xi 



Xk 



xl 



vv 



x°m 



A x 



Bs = Co 



At 


- A L 


- A L 

— A tu 




A? 


- A R 


- A R 




A L 

A st 


= Bs t 


- A L 


= Bs U 


AR 

A st 


- H R 

— D st 


- A R 


= B R 

su 



i (^+Rpo k , L 
cos 4 0w 

(flkfli ) 
(flkfli ) 

cos 2 #W 

f-R fRn"RR 

cos 2 #W 



cos' 6*w 



r 

( f L f L\2 
\J ukJ vil 



f L fL "L 

cos 2 6*w 





12 H 



cos^w 

qfe- 

•R fR\2 



3(fqkfqi 

f L f L n "L T 
R JqkJqi u ki n l 



<9 W 

f Rn"Ri 



COS 

g J qk-1 qi 



JnkJai W ki ' ft 9 



Table 8: Coefficients for x° — 



(143) 



where cy = 3 is the colour factor for gaugino decays to squarks (cy = 1 for sleptons and sneutrinos). 
And 



M xff> 



lyj (to 2 — (mj + mj') 2 )(m 2 — (m j — my) 2 ) + 4m. 2 TO 2 ,{a 2 + 6 2 } 



(144) 



where the gaugino coupling to the sfermion is parametrised as i{a(l + 75) + 6(1 — 75))}; and the 
constants a,b are for example given in |fffi| , Fig. 22-24. Note that for gaugino decays to third generation 
sfermions mixing must be taken into account by modifying the constants a,b and the sfermion mass 
appropriately. 

The gaugino decay rate to a gauge boson (W or Z) is again given by Eq. ( |143| ) , and the coupling 
constants a,b are given in [0] (Fig. 75). 

The two-body decay rates for the neutralino decays to a lighter neutralino plus a Higgs boson are 
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xi l & 


xt ud 


A s 


why 


tips? 


B s 


2(Of,) 2 


6(05) 2 


C s 






A\ 




3(/^^4i) 2 






3(/|£4i) 2 








A L 

A st 






B st 




3V2f^U kl Og 


A L 




-V2f£U kl Og 


-3V2f^U kl O[ k 



Table 9: Coefficients for x° — > X k f f 



given by 



TY-v v ffOx g 2 X(M t ,Mj,Mk)^ 
1 IX* -+ Xj^/J - 



m jk + Ff lk )(M?+Mj -Ml) 



where k=l,2,3 correspond to the three neutral Higgs h,H,A; rfu = 1 for k = 1,2 and 77I = —1 for 
fc = 3; the factors e, stand for the sign of the neutralino mass; the kinematic function A is defined by 



2j,2 



and 



ijk 



A(a, b, c)=a z + b z - c z - 4a z b 



d k 



(145) 



+ 



[MN' 2% N' 2J + M N u N[j - K^K 3 + N^)] 



2m\y sin /3 

where M,M',fi are the usual neutralino mixing parameters, and the constants c k ,d k are given by 



(146) 



ci = sin(/3 — a) , d\ = — sin(a) 
C2 = cos(/3 — a) , d 2 = cos (a) 
C3 = cos (2(3) , c?3 = cos(/3). 



(147) 



The two-body decay rates for the neutralino decays to a lighter chargino plus a charged Higgs 



boson are given by [43|: 



xfH*) 



9 m ^ M f H±Y2 [i F l + F R)(M? + Mf- Ml) + 4^^^] 



(148) 



where 



F L = cos/3K i ^i + - 7 =(^ + iV( i tan^)^2] 
F R = sin/3[^.[/ il --^=(iV^ + ^tan^)^2]. 



(149) 
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Finally the radiative two-body decays occurring through loop-diagrams have the general form: 

r «?^> = ^l^ < 150 » 

The expressions for fl^o^o can be found in Jig] . 
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C.5 Sfermion decay rates 

The most general sfermion decay widths are: 

h- \{ n ?, +l?,V m 2 - —rn 2 , — rr„-„ 



g ■ A2 [mj ,mj,mt / 

r(/i /X°) = 16 ^ m 3 — [(<& + ^4)(™| - m} - ™|o) - ±a ik b ik m f m^ k 

fi 



and 



r(/i fxj) = 16?rm 3 — y{kj + Kj){ m f. - m f- m xf} i /-//'•./"'/"" v 

where A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. 
and where: 

(a{ k \ = ( cosOf smO f ~ \ (h{ k \ (b{ k \ = ( coB0 f - sm8 f ~ \ ( f[ k \ 

h\ k = Y t (N k3 sin/3 - N M cos/?) , f Lk = sin^AT^ - ^2 (I - | sin 2 

= y * (^3 sin /3 - iV fc4 cos /?) , f t Rk = - 2 ^ sin ^(tan W iVfc2 — AT fcl ), 

/i 6 Lfc = -Y b (N k3 cos /? + N kA sin /?) , f h Lk =& sin 0^iV fcl + ^2 {\ - \ sin 2 0^) £fc , 
^Rfc = - y f> ( N k3 cos + /V fc4 sin 0) , f Rk = ^§ sin w (tan w N k2 - N kl ), 

W Lk = -Y T (N k3 cos/3 + N k4 sm(3) , f Lk = V2sm6 w N kl + V2{\- sin 2 9 W )£^, 
h Rk = ~ Y t (N k 3 cos (3 + N M sin ft) , f^ k = y/2 sin 0^ (tan 6 w N k2 - N kl ) , 
for the sfermion-fermion-neutralino interaction, and 

Zij = -Vji cos % + Y t V j2 sin % k\j = Y b U j2 cos %, 
% = Vj'i sin 8t + *t ^S'2 cos 0^ k\j = - Y b U j2 sin 9 h 

l\j = ~ U n cos °b + Y b U J2 sm 6 b , klj = Y t V j2 cos 0~ b , 
l b 2j = Uj! sin 6- b + Y b U j2 cos 9~ b , k b 2j = -Y t V j2 sin 9~ b , 

l\j = —Uji cos Of + Y T Uj 2 sin 0f , = 0, 

l 2 j = Uji sin b + Y T J7j2 cos Of, k\j = 0, 

for the sfermion-fermion-chargino interaction. 
Yf denote the Yukawa couplings: 

Yt = m t /(V2 mw sin/3), Y b = m b /(V2mw cos /3), Y r = m T /(V2 mw cos /3). 

For the first two generations the Yukawa couplings and the mixing are negligible. 
For the decay rate of the process t — ► cx° we use: 



nt L ' R ^c+x°) = 3 • 10-^(1 - -^i-^ 2 
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C.6 Two Body Decays of Sparticles to the Gravitino 



The decay width of the neutralinos to the Gravitino can be written as [21 



7 G): 

ZG) 

<PG) 



2niz T + KiZ L m h 



16vr 



k ■ m l» 
Xj_ 

16vr A| 




o 

x° 



(162) 
(163) 

(164) 



where 



Kj 7 = I cos 9\v + iV^ sin 0^ | 2 
^iZ T = |JV"a sin 9w — N' i2 cos ^w\ 2 
Kiz L = | A^ 3 cos P ~ N 'a sin /3| 2 
K ih o = \N- 3 sin a — A^ 4 cos a| 2 
Kj^o = | JV^3 cos a + N' iA sin a| 2 
«M« = \ N k sin /? + cos /3| 2 



and </? = (h°, H°, A ) is any of the neutral Higgs scalar bosons. The couplings N-j are given in 
The chargino decay widths into gravitino final states are given by 



(165) 
Eq. ©. 



VF+G) 

# + G) : 



in 



w 



16vr 

K iH+ xj 
16vr A| 



A 4 



xj 



m 



H+ 



m. 



(166) 
(167) 



with 



KiW L 
K iH+ 



5 0*1 



+ 1*7, 



(1 



|y i2 | 2 sin 2 /3+|[/ j2 | 2 cos 2 /? 
|^ 2 | 2 cos 2 /3 + |C/ i2 | 2 sin 2 /3 



and the widths of the 2-body slepton decays into gravitinos are given by 

r(/-/G) 



i 



8vrA 4 



(168) 



(169) 



C.7 R p violating two body decays 



For the Rp violating Yukawa coupling \'l- k the decay rate of a squark q~\ of flavour i into the final 
states qu — ► qjgfc is given by 



li 



A" j 2 sin 2 0g|pj 



8ttM 



2fe -Pfc) 



(170) 
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where 9 q is the mixing angle of the squarks qu, qm which form the mass eigenstates qu, qu. Here 

\Pi\ = ((ML ~ (M qj + M qk f)(M? u - (M gj - M qk ) 2 ))l/(2M qii ) 
( PrPk ) = ((Ml + \p i \ 2 )(M 2 qj + \p l \ 2 )) 1 2+\p l \ 2 (171) 

The decay rate for q2i — ¥ Qj<lk can be simply obtained from the above formula by replacing sin#q —* 
cos6g and qu — > q2i- The decay rates for the other operators are analogous, except that the mixing 
term sin 9 has to be replaced by cos 6 for those sfermions which are part of the SU(2) doublet superfields 
(the Li j and Qj fields in Eq. (|l|)). Note also that the result has to be multiplied by a colour factor 
Cf = 3 for the slepton decays via LiQjD k . 

C.8 R p violating three body decays 

Neutralinos can decay through R p violating couplings to SM particles via Eq. (p9[)-(j3l|). The Matrix 
element squared for the decay x^i ~~ > SiUjd k of the operator LiQjD k is given by pq| : 



\M($ -+ e iUj d k )\ 2 =8c f g 2 y 2 k { (172) 
D(uj) 2 ei ■ d k [(a(uj) 2 + b(u) 2 )xj ■ Uj + 2a( Uj )b(u)m Uj M^o ] 
+ D(d k ) 2 ei ■ Uj [(a(d k f + b(d) 2 )^ ■ d k - 2a(d k )b(d)m dk M^ } 
+ D(ej) 2 Uj ■ d k [(a(ei) 2 + 6(e) 2 )x° ■ (k + 2a(ei)b(e)m e M^ ] 

- D(ei)D(uj) [a(uj)a(ei)m ei m Uj Xi ■ d k + a(u j )b(e)in Uj M jt oe i ■ d k 

+a(ei)b(u)m ei M i0 Uj ■ d k + b(e)b(u) Q(uj ,Xi , (k,d k )] 

- D(uj)D(d k ) [a(u j )a(d k )m Uj m dk x < i • &i - a(uj)b(d)m Uj M^oei ■ d k 

+a(d k )b(u)m dk M^oei ■ uj - b(u)b(d)Q(uj,Xi,d k ,ei)] 

- D(ei)D(d k ) [-a(ei)b(d)m et M^oUj ■ d k + a(ei)a(d fc )m ei m rffe x[ ) ■ Uj 

+a(d k )b(e)m dk M jt oe i ■ uj - b(e)b(d) G(Xi, <H, Uj, d k )]\ 

where the function Q(a, b, c, d) = (a-b)(c - d) — (a-c)(b-d) + (a- d)(b- c). Here cj = 3 is the colour factor 
and g is the weak coupling constant g = . \ ■ We have denoted the 4- momenta of the initial and 
final state particles by their particle symbols. The functions D(pi) denote the propagators squared 
for particle pi and are given by 

D(p i )- 1 = M| ? + w 2 . - 2x d ■ Pi - m 2 pr (173) 

The coupling constants a(pi), b(p) are given in Table [U]. 

The amplitude squared of the decay to the neutrino, x? — * v idjd k , can be obtained from the 
above result by a set of simple transformations of the 4-momenta, the propagator functions D and the 
couplings a(pi), b(p) : — > z/j, Uj — > dj. 

The result for the operators LiLjE k is completely analogous, except the colour factor Cf = 1. The 
matrix element squared for the operator UiDjD k is 

\M(xi - n44)| 2 =8c f g 2 X" 2 { (174) 
D(ui) 2 dj ■ d k [(b(u) 2 + a(ui) 2 )x° ■ «i - 2b(u)a(ui)m Ui M^] 
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Ui 

di 



TV 7- " 

JV 3i 



2M W cos /3 m ei 



2M W cos /3 me » 




N' 

_4i 



2Mi V sin/3 m «j 



2Mi V sin/3 m «j 



A': 



31_ 



2Mvk cos (3 md k 



2M W cos /3 md k 



b(Pi) 
- tan e w N' u 



\(N' 2l -tan^^i 



1/7 



■i(^-itau^JVj { ) 



-± tan w i\r{, 



Table 10: The Coupling constants a(pi) and used in the LSP decay calculation for a neutralino 



xl 



+D(d j ) 2 u i ■ d k [(a(djY + 6(d) )Xi • 4 " 26(d)a(d i WAfgo 



+D(d k ) 2 Ui ■ dj 
+D{ui)D{d k ) 



[(a(4) 2 + &(d~) 2 )x° • 4 - 2b(d)a{d k )m dk M 9 o 



j Xi- 



+D{ui)D{d k ) 
+D(d J )D(d k ) 



[a(ui)a(dj)G{ui,x°i, d ji d k) + K u )K d ) m u i m d k Xl ■ 4 
-b(u)a(dj)m Ui M^odj ■ d k - b{d)a{ui)m dj M^}Ui ■ d k ] 

[a(ui)a(d k )g(ui,xi,d k ,dj) + b(u)b(d)m Ui m dk xi • 4 
-b(u)a(d k )m Ui M^,dj ■ d k - b(d)a{ui)m dk M^oUi ■ dj] 

[a(dj)a(d k )g(dj,Xi,d k ,Ui) + b{d)b(d)m dj m dk x® - Ui 

-b(d)a(d k )m dj M^oUi • J fc -b(d)a(dj)m dk M^oUi ■ dj]\ 

Here the colour factor c/ = 6. The partial decay rates for the neutralino decay modes is then given by 

T^{LLE) = J ^ T ^_(2|^(x[ ) - e % vfe k )\ 2 + 2\M(x° -> u^e^dE^ 

r^o (LQD) = J^ T ^(2|A^(x[ ) ^e l n i 4)| 2 + 2|Ai(x[ ) ^M J 4)| 2 )^^ i (175) 

T^UDD) = S^^^MijQ^Uidjd^dEidEi 



where the matrix element squares are multiplied by a factor of two since the LSP is a Majorana 
fermion and can decay to the conjugate final states. 



Charginos can decay via Eq. (32)-(|3j), and the matrix elements squared for the LiQjD k operator 

(176) 



can be written as [22]: 
\M(xf mjdk)\* 



4c /ff 2 A' 2 



a. 



R 2 (e iL ) 



iXi ■ "i)( u j ■ d k ) 



+ 



(n ■ d k ) 
R 2 (d jL ) 



'1 + PR)(xt ■ Uj) + 2Ke(l3 L {3* R m uj M~ t )} 
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-He I 



{R(e iL )R(d jL ) v Ai 



^ Mid*)! 5 



4c /5 2 A' 2 



(dj ■ d k ) 



(e% ■ d k ) 



#j7^r {ill + llMt ■ ei) + 2Ke( lLl * R m ei M~ t )} 

~ R{j> L )R{u L ) ne { lh5 * L ^^i' e *' ^ + lL5 R m dj M~+(ei ■ d k ) 
+-/RS* L m ei M^+(dj • 4) + lRS R m e im dj (xi' ■ 4)}] 



/2 1 1 



\M{xf eiUjUfc)] 2 
|A<(Xi + Mj« fc )| 2 



cos* PR*(d kR ) { jAXl k) 



(177) 
(178) 

(179) 



where the couplings a L ,R, (3l,r,7l,r and 5 l ,r are 



a L = 0, 
— ' 

7L = ^, 



a R = -iUn 

= igm ei U l2 
,H V2M W cos /3 

X „, X igm d jUi2 



The final state momenta are denoted by the particle symbols. is the charg ino mass and i^iei,dj,dk 

are the final state fermion masses, c/ = 3 is the colour factor. The propagator terms R(p) are defined 
by 



R{e lL ) 


= (x + 




~ ™>liLi 


R(d jL ) 


= (x + - u,r - 


R&l) 


= (x + 


-e t f 


~ ™>liL> 


R(u jL ) 


= (X + - d 3 ) 2 - 


R(dk R ) 


= (x + 


-u k ) 2 


~ 2 
~ m dkRi 


R(e kR ) 


= (X + -Vkf- 


R(UiR) 


= (x + 


- dtf 


~ m uiRi 


R(d jR ) 


= {x + ~u j f- 



™2 



™lkR 
™2 



For the operator LiLjE k the matrix elements squared are 



\M{xt -► eie jek f 



4g 2 X 2 



ixt -Vj)iyi -e k ) 



Tie 

.2 \2 



^-—g(p,Vi,e k ,Vj 



(180) 



-4 



R(e iL )R(e jL )' 

{ill + t|)(x/ • ei) + 2fte( 7L7 £m ei M- j+ )} 

^T^y + • ej) + 27*e(* L ^m ci M-+)} 

~ R{v L)R{B L ) ne {^^(xt^^k, ej) + -y L 6 R m ej M~+( ei ■ e k ) 
+jR5* L m ei M^+(e j ■ e k ) + ^ R S R m ei m ej (xi' ■ e fc )| (181) 
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\M(xt ^-i/fc)| 2 



Va 2 


4 


_^ 2 (e fc i?) 


Va 2 


4 


i? 2 (efcii) 



(182) 
(183) 



where a, (3, 7, (5 are given as above except that in 5r m<« is replaced by m e j and = because of 
vanishing neutrino mass, and = — J^^ ek ^ 2 s/3 - ^ or ^ e P era ^ or UDD the matrix elements squared 
are 



\M(xt ^ did 3 d k )\ 2 = 4c f g 2 X 2 
\M{xt u iUj d k )\ 2 = Ac f g 2 \ 2 



\M(xt -> Uid jUk )\ 2 = 4c f g 2 X 2 



ft 



(X + • di){dj ■ d k ) 



UiR) 



R 2 (d 



• Uj)(ui ■ d k ) 



R 2 (d 



■jR) 
5 2 

— (x + ■ t*fc)(t*i ■ d,0 



Pr,1r,Sr are given by 



PR = * 



Sr 



V2M w sin(3 
igm dk U l2 



_ igm dj Ui2 
1R V2M W cos p 



\f2M w cos/3 

and the colour factor cf = 6. The partial decay rate for these modes is then given by 



T-+(LLE) 



T~+(LQD) 



r -+(c/£»l>) 



/ 2^ WT( I^(Xj + Wfc)| 2 + |A4(x+ eie jefe )| 2 + 

|M(x/ e^^fc)! 2 + l-A^CXi 1- ViejVk^dEidEj 

|7W(X Z + -► ej^Ufe)! 2 + |.M(x/ -► VidjU k )\ 2 )dEidEj 
( |M(x+ - J4-4))| 2 + - «i«jd fc )| 2 + 



r 1 L_ 



\M{xt UidjU^dEidEj 



(184) 
(185) 
(186) 



(187) 



(188) 
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D Common Blocks 



In the following the common blocks are listed in alphabetical order. 

• COMMON/ CHAMIX/ OIJL(2,2), OIJR(2,2), V(2,2), U(2,2), FM(2), ETA(2) 

• COMMON/ CHANEU/ OIJLP(4,2), OIJRP(4,2) 
OIJL, OIJR : chargino couplings to Z. 

V, U : diagonalising matrices of chargino.. 
FM : positive chargino masses. 
ETA : sign of diagonal values. 

OIJLP, OIJRP : chargino couplings to neutralinos. 

• COMMON/ CONST/ IDBG, IGENER, IRAD, FRAD: (Input cards DEBUG, GENER, ISR and 
FSR) 

• COMMON/ COUPLS/ A,B: Coupling constants used in the R p violating gaugino decay calcu- 
lation. 

• COMMON/ DECSEL/ IDECSEL(23),IPROCSEL(ll): (Input cards DECSEL,PROCSEL) 

• COMMON/ DOTPROD/PP: Dot-products used in the R p violating gaugino decay calculation. 

• COMMON/ FINDEX/ FMPR1, FMPR2, XCROST, APRO 
FMPR1, FMPR2 : masses of the 2 current produced sparticles. 
XCROST : current cross section. 

APRO : maximum differential amplitude. 

• COMMON/ GAUG2BDECS/ IJKGAUG(MAXGAUG2BODDECS,3),NGAUG 

• COMMON/ GAUG2BDECS2 / GIJKGAUG(MAXGAUG2BODDECS) 
NGAUG : Number of R p conserving 2-body gaugino decay modes. 

IJKGAUG(n,i) : i=1...3 gives the initial(i=l) and final state particles (i=2-3) in the 2-body 
decay. 

GIJKGAUG : Width of 2-particle decay. 

• COMMON/ GAUGM/ MASSNEU,MASSCHA : Neutralino and chargino mass used in current 
R p violating gaugino calculation. 

• COMMON/ HABKANE/XM(i,j) : Mixing matrix N' tj defined by Eq. ©. 

• COMMON/ HIGGSES/MLH,MA,MHH,MHPM,SINA,COSA : Higgs masses and mixing angles. 

• COMMON/ INDEXX/ INDEX, INDEX1, INDEX2, NEVT 

INDEX is an internal variable determining whether the program currently processes neutralino 
(INDEX=1), chargino (INDEX=2) sfermion (INDEX=3), gravitino (INDEX=4), single spar- 
ticle production (INDEX=5), or Higgs (INDEX=5). 

INDEX1, INDEX2 are the index code of the specific type of 2 sparticles produced. 

• COMMON/ INDICES/I,IJ,IJC,IJF,IJF_MIX 
1(4) : Input Card INDIC 

IJ(i=1..3,j=1..4) : Decay product(i) of R p violating neutralino(j) decay. 
IJC(i=1..3,j=1..2) : Decay product(i) of R p violating chargino(j) decay. 

IJF(n=1..68,c=1..2,fs=1..2) : Decay products of direct R p violating sfermion decays. n=left- 
handed sfermion decaying into the final state fs=1..2. There are up to c=1..2 decay modes. 
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IJF_MIX : Same as IJF, except that the n=right-handed sfermion. 

COMMON/ ISR/ QK(4) (4-vector of ISR 7). 

COMMON/ KINEM/ FLUM, ECM, S, ROOTS, T, Q, Q2, EN(2) 

FLUM, ECM, S, ROOTS : L, E C m, s, yfs ( L is the luminosity, E CM is the nominal energy, 

s and yfs change when initial state radiation is permitted). 
T, Q, Q2, EN : t,q,q 2 and energies of 2 products. 
COMMON/ LAMDA/XLAMA : R p violating coupling strength A. 
COMMON/ LIFE/LIFETIME : Input Card LIFE. 

COMMON/ LSUSY/XLSUSY : Spontaneous Supersymmetry breaking scale Xsusy- Used to 
calculate the mass of the gravitino. 

COMMON/ LUSHOWERS/LISTQMX 

COMMON/ LUSHOWERS2/LISTSHOWERIND,NSHOWERS 

LISTQMX(i),LISTSHOWERIND(i,j) : Give the invariant mass and the index of the par- 
ticles (j) of the i-th shower which are to be processed by JETSET routine LUSHOW. 
NSHOWERS : Number of showers to process. 
COMMON/ MDS/MODES,MIX : Input Cards MODES and MIX. 
COMMON/ MSSM/ TANB, SINB, COSB, FMGAUG, FMR, FMO, ATRI(3) 
TANB, SINB, COSB : tan (3 = v 2 /v 1. 
FMGAUG : M 2 gaugino mass. 
FMR : Higgs mixing term \i. 
FMO : common scalar mass txiq. 

ATRI : 3 trilinear couplings for stop, sbottom and stau respectively. 

COMMON/ NEUMIX/ ZR(4,4), WAS(4), ESA(4), VOIJL(4,4), VOIJR(4,4), GFIR(4,4), GFIL(4,4) 

ZR : neutralino diagonalising matrix. 

ESA : phase factors of eigenvalues . 

WAS : absolute masses in ascending order. 

VOIJL, VOIJR : neutralino couplings to Z. 

GFIL, GFIR : neutralino couplings to u, d, i>, e. 

COMMON/ REORDER/ ISPA(12,2), KL(2,18), KLAP(2,18), IDECS(12,2) (Correspondence 
matrices for decays and equivalence with LUND) 

ISPA : LUND codes of left/right sfermions (18 possible). 

KL : standard particles accompanying a gaugino in 3-body decay . 

KLAP : IDs of t,u exchanged particles. 

IDECS : mapping of SUSYGEN standard particle convention to the LUND ones. 

COMMON/ RESA/RES,RES1A,RES2A : 

RES : Total Decay rate of R p violating neutralino decay. 

RES1A,RES2A : Partial Decay widths of R p violating neutralino decays to final states of type 
(1),(2). 

COMMON/ RKEY/rgmaum, rgmaur, rgmO, rgtanb, rgatri, rgma, rfmsq, rfmstopl, rfmstopr, 
rfmsell, rfmselr, rfmsnu, rfmglu, recm, rfTum 

rgmaum,..,rfmglu : The single precision version of the variables in the COMMON STEER. 
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recm,rflum,rvscan : CM-Energy, Luminosity in single precision. 

COMMON/ RPARI/RPAR,MGM 

RPAR : Logical, TRUE — ► R p is conserved. 

MGM : Logical, TRUE — > Gauge Mediated SUSY gravitino decays enabled. 

COMMON/ RWIDTHS/GIJ,GIJC,GIJF,GIJF_MIX : Decay rates of the R p violating gaugino 
and sfermion decays. See also COMMON INDICES. 

COMMON/ SDECAY/ DAS, DBS, DCS, DATL, DAUL, DATUL, DASTL, DBSTL, DASUL, 
DBSUL, DATR, DAUR, DATUR, DASTR, DBSTR, DASUR, DBSUR, XDEC (17,64) 

DIMENSION CURENT(17) 

EQUIVALENCE (CURENT(1),DAS) 

DAS etc correspond to the 17 variables that characterise gaugino 3-body decay. There are 64 
different patterns; 

1. 6x4 neutralino-neutralino to uu, dd, vv, 11. 

2. 8x4 chargino-neutralino to ud, lv. 

3. 1x8 chargino-chargino uu,dd, vv, 11. 

COMMON/ SM/ FMW, FMZ, GAMMAZ, SINW, COSW, ALPHA, E2, G2, PI, TWOPI, 
FLC(12), FRC(12), GMS(12), ECHAR(12) 

FMW : Wmass (80.2GeV). 

FMZ : Z mass (91.19GeV). 

GAMMAZ : Z width (2.497GeV). 

SINW, COSW : sin W (0.231243), cos/% Weinberg. 

PI, ALPHA, e2, G2 : vr, ±, e 2 = 4vra, g 2 = 

FLC : T 3 - Qsin 2 6 W , left SM couplings. 

FRC : -Qsm 2 6 w , right SM couplings. 

GMS : masses of standard particles in the following order (up, down, neutrino, electron) x3 for 

the three families. 
ECHAR : charge of standard particles. 

COMMON/ SPARC/ ZINO, WINO, SELE, SMUO, STAU, SNU, SQUA, STOPA, SBOTA (In- 
put cards ZINO etc) 

COMMON/ SPARTCL/ FMAL(12), FMAR(12), RATQA(12), FGAMC(12), FGAMCR(12), 
COSMI(12), FM1(12), FM2(12) 

FMAL(12), FMAR(12) : Left/right masses of sparticles. 
RATQA(12) : charge. 

FGAMC(12), FGAMCR(12) : left/right coupling to Z. 
COSMI(12) : cosine of mixing angle (cos</>). 

FM1(12), FM2(12) : masses of the 1,2 combinations for the 12 sfermions. 
All the above are in the following order: u,d,u,e x3 generations. 
COMMON/ SSCALE/RSCALE : Input Cards RS. 

COMMON/ SSMODE/ NSSMOD, ISSMOD(MXSS), JSSMOD(5,MXSS) 

COMMON/ SSMOD1/ GSSMOD(MXSS), BSSMOD(MXSS) (ISAJET type commons, for com- 
parison) 
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MXSS : maximum number of modes. 
NSSMOD : number of modes. 
ISSMOD : initial particle. 
JSSMOD : final particles. 
GSSMOD : width. 
BSSMOD : branching ratio. 

COMMON/ STEER/ GMAUM, GMAUR, GMO, GTANB, GATRI, FMSQ, FMSTOPL, FM- 
STOPR, FMSELL, FMSELR, FMSNU, FMGLU 

(M, fi, mo, tan/3, A, m^, , , m e ~ L , m^, m^, nig) 

COMMON/ STOPMIX/STOPl,SBOTl,STAUl,STOP2,SBOT2,STAU2, PHIMIX1,PHIMIX2,PHIMIX3 
: Input cards MASMX 

COMMON/ STR/ WRT, SCAN, LEPI (Input cards LUWRIT, SCAN, LEPI) 

COMMON/ SUSY/TANTHW,COSB,SINB : tan 9 W , cos 0, sin 0, where is the ratio of the neu- 
tral Higgs vevs. 

COMMON/ TWOBODIES/ TWOB : Input cards TWOB 
COMMON/ UBRA/ NDECA(-80:80) 

NDECA : number of decay nodes open to the particle with corresponding ID. 

COMMON/ UBRA1 / BRTOT(2, 50, -80:80) 

BRTOT : i)pointer to SSMODE, ii)integrated branching ratios. 

COMMON/ WIDTHS/ GW, WIDFL(12), WIDFR(12), FMELTW, FMERTW, FMELUW, FMERUW, 
LINDA(18,6,6) 

GW : The product of Z (or W) mass and width (GeV 2 ). 
WIDFL, WIDFR : Sparticle widths (GeV). 

FMELTW /FMERTW, FMELUW /FMERUW : The products of sparticle masses and 

widths {GeV 2 ). 
LINDA : (Internal use). 

COMMON/ VARIABLES/ FMS, FMI, FMK, FML1, FML2, ETAI, ETAK, BRSPA(6,48), LIND(6,6,6), 
BRGAUG(23,6,6), FMELT, FMERT, FMELU, FMERU 

All variables relate to the R-parity conserving gaugino decays except for BRSPA, which relates 
to the R-parity conserving sfermion decays. 

FMS : is the s channel exchanged particle W or Z. 
FMI, ETAI : is the mass,7? of the father. 
FMK, ETAK, : is the mass,r/ of the son. 
FML1, FML2 : are the accompanying partons. 

FMELT/FELU, FMERT /FMERU : are the sparticles exchanged in u,t channels. 

— LIND(i, j, k) are the permitted transitions; 

i : is one of the 6 patterns uu, dd, 11, vv, ud, lv. 
j : is one of the son gauginos. 
k : is one of the father gauginos. 

- BRGAUG(i, j, k) are the integrated differential widths; 

i : is one of the 18 patterns (uu, dd, vv, 11, ud, lv) x3 generations. 
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j : is one of the son gauginos. 
k : is one of the father gauginos. 
— BRSPA(ij) are the partial sparticle widths; 
i : is the index of the daugther gaugino. 
j : is the index of the father sfermion. 

• COMMON/ XCROS/ XGAUG(8), XETA(8) 

XGAUG are the masses of the 4 neutralinos and the 2 charginos %i 
XETA the rj = ±1 phase factors. 
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Ntuple variables 



Variable 


Description 


Variable 


Description 


Variable 


Description 


fECM 




fM2 


M 2 


fmu 


A 4 


fmO 


m 


ftanb 


tan/3 


fA 


A t 


fmA 


m A 











Table 11: Input parameters - block "COL1". 



Variable 


Description 


Variable 


Description 


Variable 


Description 


fxOl 


M(xS) 


fx02 


M(x u 2 ) 


fx03 


M( X U 3) 


fx04 


M(xl) 


fxlp 


M(xt) 


fx2p 


M(xt) 


fupl 


M(u L ) 


fdownl 


M(d L ) 


fneu 


M(y) 


fel 


M{e L ) 


fupr 


M{u R ) 


fdownr 


M(d L ) 


fer 


M{e R ) 


ftopl 


M{h) 


fbotl 


M(6i) 


ftaul 


M(n) 


ftop2 


M(t 2 ) 


fbot2 


M(b 2 ) 


ftau2 


M(f 2 ) 


fmhl 


M(h) 


fmh2 


M(H) 


fmhpc 













Table 12: Masses - block "COL1". 



Variable 


Description 


Variable 


Description 


Variable 


Description 


zxOla 




zxOlb 


N 2 i 


zxOlc 


iVsi 


zxOld 


N 41 


zx02a 


N 12 


zx02b 


^22 


zx02c 


A^32 


zx02d 


N 42 


zx03a 


#13 


zx03b 


^23 


zx03c 


N33 


zx03d 


N 43 


zx04a 


N U 


zx04b 


N 2A 


zx04c 


N u 


zx04d 


N u 


ull 


Un 


ul2 


Ui 2 


vll 


Vn 


vl2 


Vn 







Table 13: Neutralino/Chargino eigenvectors - block "COL3". 
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Variable 


Description 


Variable 


Description 


Variable 


Description 


sxlxl 




sxlx2 




sx2x2 




sxlx3 


'(x?X§) 


sx2x3 


<xlxl) 


sx3x3 




sxlx4 


'CM 


sx2x4 




sx3x4 


*(X°3Xi) 


sx4x4 


<xlxl) 


sxlxlp 




sxlx2p 


a (xix 2 ) 


sx2x2p 


v{xixi) 


snue 


o{y e v e ) 


snutot 


a(u e u e ) + oiy^Vp) 
+ a(v T v T ) 


ssel 


cr(e L e L ) 


sser 


v{e R e R ) 


sselser 


cr(e L e R ) 


smul 




smur 




staul 


°"(fifi) 


stau2 


0"(f 2 f 2 ) 


sbtl 


<t(&i&i) 


sbt2 


a(b 2 b 2 ) 


stpl 




stp2 


v(t 2 t 2 ) 


sqlr 




sgxl 


<X\G) 


sgx2 


v(x° 2 G) 


sgx3 


<xl&) 


sgx4 


°{xlG) 


srpvxl 


<X i») 


srpvx2 




srpvx3 


<X» 


srpvx4 


<X» 


srpvxlp 




srpvx2p 


*(xtn 

a(hA) 


szhl 


a{hZ) 


szh2 


a{HZ) 


sahl 


sah2 


a(HA) 


shphp 


a(H + H~) 



Table 14: Cross sections in pb - block "COL2". 



Variable 


Description 


Variable 


Description 


Variable 


Description 


sxOl 


I\,(Xi) 


sx02 


?rM) 


sx03 


r Rp (xi) 


sx04 


i\,(x2) 


sxpl 


FrM) 


sxp2 


r Rp (xi) 


sul 


t R p (ul) 


sell 


T Rp (d L ) 


snel 




sel 


FR p (e L ) 


scl 


Fr p {c l ) 


ssl 


Tr p (sl) 


snmcl 


T R P (^) 


sml 


t r p (^l) 


stl 




sbl 




sntl 


T R P {Vj) 


stal 


r Rp (fi) 


sur 




sdr 


TrMr) 


sner 




ser 


t R p (.zr) 


scr 


t R p (cr) 


ssr 


t R p (sr) 


snmer 




smr 




str 


r Rp (t 2 ) 


sbr 


T Rp {b 2 ) 


sntr 




star 


r i? P (^2 



Table 15: R-parity conserving widths (Tr p ) of sparticles in GeV - block "COL4". Does not include 
GMSB decays to the gravitino. 
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Variable 


Description 


Variable 


Description 


Variable 


Description 


rxOl 


r%,(xi) 


rx02 


r&,(x 2 ) 


rx03 




rx04 




rxpl 




rxp2 




rul 




rdl 




rnel 




rel 




rcl 


%(cl) 


rsl 




rnmel 




rml 




rtl 


rbl 


%(&i) 


rntl 




rtal 


%(n) 


rur 




rdr 




rner 




rer 




rcr 




rsr 




rnmer 




rmr 




rtr 


%(i 2 ) 


rbr 


%(6 2 ) 


rntr 




rtar 





Table 16: R-parity violating widths of sparticles in GeV - block "COL5". 



Variable 


Description 


Variable 


Description 


Variable 


Description 


rxOl 


^gmsb(xi) 


rx02 


^gmsb{xV) 


rx03 


^gmsb{xV) 


rx04 


Fgmsb(xI) 


rxpl 


^GMSBixt) 


rxp2 


^GMSsixt) 


rul 


^gmsb(ul) 


rdl 


^GMSB{dl) 


rnel 


^GMSB^e) 


rel 


^gmsb(^l) 


rcl 


Tgmsb(cl) 


rsl 


Fgmsb(sl) 


rnmel 


^GMSb{^h) 


rml 


FgmsbOj-l) 


rtl 


^GMSB(h) 


rbl 


^GMSB{b\) 


rntl 


^gmsb{vj) 


rtal 


Fgmsb(ti) 


rur 


^gmsb{ur) 


rdr 


^GMSBidn) 


rner 




rer 


^gmsb(gr) 


rcr 


^gmsb{cr) 


rsr 


Fgmsb(sr) 


rnmer 




rmr 




rtr 


^GMSB(h) 


rbr 


^GMSBifa) 


rntr 




rtar 


Tgmsb{t2) 


Table 17: GMSB widths (Tr gmsb ) 


of sparticles in GeV - block "COL6". 


Variable 


Description 


Variable 


Description 


Variable 


Description 


rxOl 


r(x!f) 


rx02 


r(xS) 


rx03 


r(x8) 


rx04 


r(x8) 


rxpl 


nxt) 


rxp2 


r(x 2 + ) 


rul 


r(«L) 


rdl 


nd L ) 


rnel 


m) 


rel 


r(e L ) 


rcl 


r(s L ) 


rsl 




rnmel 


r(^) 


rml 


r(£L) 


rtl 


r(ti) 


rbl 


r(5i) 


rntl 


r(* T ) 


rtal 


r(n) 


rur 


r(«ji) 


rdr 




rner 




rer 


r(e«) 


rcr 




rsr 




rnmer 




rmr 




rtr 




rbr 


r(6 2 ) 


rntr 




rtar 


r(f 2 ) 



Table 18: Total widths (r = r Rp + % + r RGMgB ) of sparticles in GeV - block "COL7". This block 
is only present in the GMSB or /? p modes. 
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Variable 


Descript 


ion 


Variable 


Description 


brupxOl 


dT(u L — 


> ^Xi) 


brupx02 


dT(u L — 


' ^X2) 


brupxpl 


dT(uL — 


^xf) 


brdownxOl 


dT{d L - 


dx?) 


brdownx02 


dT{d L - 




brdownxpl 


dT{d L - 


u Xi ) 


brneuxOl 


dT{v e 


^eX?) 


brneux02 


dT(v e -> 


^eX?) 


brneuxpl 


al [y e — > 


eT Xi ) 


brelxOl 


dT(e L -» 


ex?) 


brelx02 


dT(e L - 




brelxpl 


dr(e L -» 


^eX^) 


brstxOl 


dT{ti -► 


c(ort)xi) 


brstx02 


dT(ti -► 


c(ort)x?) 


brstxpl 


dr lh -► 


&xr) 


brbotxOl 


dT(&i -► 


fox?) 


brbotx02 






brbotxpl 


dT(6i -»• 


c{ort)x 2 ) 


brntxOl 


dr(p T -> 


^X?) 


brntx02 


cir(z> T — ► 




brntxpl 


dr(p r — > 


r^) 


brtauxOl 


dT(fi - 


rx?) 


brtaux02 


dT(fi - 




brtauxpl 


dT(fi 


^rX^) 



Table 19: R-parity conserving partial widths in GeV - left handed sfermions. Block "COL8". 



Variable 


Description 


Variable 


Description 


brupxOlr 

brupxpl r 

brdownx02r 

brelxOlr 

brelxpl r 

brstx02r 

brbotxOlr 

brbotxplr 

brtaux02r 


dT(u R -> u\i) 
dT(u R -> dxf) 
dT(d R - dx^) 
dr(e R - ex?) 
dr(e R -> f e Xi~) 
dT(* 2 c(ort)x?) 
dr(6 2 - 6x?) 
dT(b 2 ^ c(ort) X2 ) 
dF(f 2 -+ rx?) 


brupx02r 

brdownxOl r 

brdownxpl r 

brelx02r 

brstxOlr 

brstxplr 

brbotx02r 

brtauxOlr 

brtauxplr 


dT(u R -> ux 2 ) 

dr(d fl - d x ?) 

dr(J R -> nxr) 
dT(e R - e X ^) 
dr(t 2 - c(ort)x?) 

dr(i 2 - & x +) 
dr(6 2 - & x °) 

dr(f 2 - rx?) 

dr(f 2 - ^) 



Table 20: R-parity conserving partial widths in GeV - right handed sfermions. Block "COL8". 



Variable 


Description 


Variable 


Description 


brx2a 


dT( X u 2 - 


QQXi) 


brx2b 


dTixt - 


e+e xi) 


brx2c 




M + M~X?) 


brx2d 


dT(x° 2 - 


r + r-x?) 


brx2e 


dr( x °2 - 


vvxi) 


brx2f 


dT(x° 2 - 


7X?) 


brx2g 


dT(x° 2 - 


h(orH)x<() 


brx2h 


dr(x§ - 


H^xf) 


brxpla 


dr(xt - 


qq'xi) 


brxplb 


dr(x^ - 


evxi) 


brxplc 


dT(xt - 


V-vXi) 


brxpld 


dr(x^ - 


r^X?) 


brxple 


dT(xf ~ 


H+ X \) 


brxplf 






brxplg 






brxplh 







Table 21: R-parity conserving partial widths in GeV - Charginos and Neutralinos. Block "COL9". 
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Variable 


Description 


Variable 


Description 


Variable 


Description 


exlxl 
ex 1x3 
exlx4 
ex4x4 
ex2x2p 

esel 

emul 

etau2 

etpl 

edowl 

eupr 

estr 

egx3 

erpvx2 

erpvxlp 

ezh2 

ehphp 


e(xlx°3) 
eixlxl) 
eUlxl) 

<xixt) 

e(T 2 f 2 ) 

eiiih) 

e(d L d L ) 

e(u R u R ) 

e(sRS R ) 

<X° 3 G) 

e(x°2") 

<xtn 

e(HZ) 
e{H+H-) 


exlx2 

ex2x3 

ex2x4 

exlxlp 

enue 

eser 

emur 

ebtl 

etp2 

echal 

edowr 

egxl 

egx4 

erpvx3 

erpvx2p 

eahl 


<xlxl) 
<X°2XI) 

4xtxt) 

e(z> e P e ) 
e(e R e R ) 

e{d R d R ) 

<X\G) 

e{xlG) 

<X» 

e(xP~) 
e{hA) 


ex2x2 

ex3x3 

ex3x4 

exlx2p 

enutot 

eselser 

etaul 

ebt2 

eupl 

estl 

echar 

egx2 

erpvxl 

erpvx4 

ezhl 

eah2 


<X§X§) 

<xlxl) 
tixtxt) 

e(u e u e + 

+v T v T ) 

t(e L e R ) 

e(fifi) 

e(b 2 b 2 ) 

t(s L s L ) 
t{c R c R ) 
<X%G) 

e(hZ) 
e(HA) 



Table 22: Efficiencies (calculated in routine USER) for the various signals - block "COL10". 
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F Test Run Input 



**** SUGRA mode 

MODES 1 

M 80. 

mu -50 . 

mO 100. 

tanb 1.41 

ECM 183. 

* initial state radiation ON 
ISR 1 

* list events 
DEBUG 1 

* write o/p to fortran file unit 12 
GENER 1 

LUWRIT TRUE 

**** generate one chargino pair event 
WINO TRUE 
SUSEVENTS 1. 
END 



G Test Run Output 
G.l The susygen.dat file 



INPUTS : 

M 

mO 

At 

Ecm 



80 . 000 mu 

100.000 TANB 

. 000 Ab 

183.000 EVENTS 



Sparticle masses 



SUPR 

SDNR 

SELR 

SNU 

STPL 

SBTL 

STAL 



253 . SUPL 

253. SDNL 

110. SELL 
117. 

312. STPR 

263 . SBTR 

125. STAR 



260. 
263. 
125. 

307. 
253. 
110. 



-50.000 
1.410 

0.000 Atau = 0.000 
1.000 RAD C0RR= 1 



Ml = 40.107 M2 = 80.000 

NEUTRALINO m, CP, ph/zi/ha/hb 1 = 

NEUTRALINO m, CP, ph/zi/ha/hb 2 = 

NEUTRALINO m, CP, ph/zi/ha/hb 3 = 

NEUTRALINO m, CP, ph/zi/ha/hb 4 = 

CHARGINO MASSES = 82.205 

CHARGINO ETA = 1.000 



45.0 1. 0.755 -0.196 -0.083 -0.620 

49.7 1. 0.625 0.014 -0.118 0.772 

98.9-1. 0.054 -0.476 0.872 0.099 

124.4 1. 0.192 0.857 0.467 -0.100 

122.500 
1.000 
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U matrix WINO HIGGSINO 
W1SS+ -0.076 0.997 

W2SS+ 0.997 0.076 



V matrix WINO HIGGSINO 
W1SS- 0.722 -0.692 

W2SS- 0.692 0.722 



HIGGS masses 



Light CP-even Higgs = 56.372 

Heavy CP-even Higgs = 313.893 

CP-odd Higgs = 300.000 

Charged Higgs = 310.070 

sin(a-b) = -0.613 

cos(a-b) = 0.790 

PARENT — > DAUGHTERS WIDTH (eV) BRANCHING RATIO 



MASS 


of Z1SS 




44.9614410400391 


GeV 






Z2SS 


— > 


Z1SS 


UP 


UB 


0.019 


2. 


, 12 


Z2SS 


— > 


Z1SS 


DN 


DB 


0.024 


2 


.73 


Z2SS 


— > 


Z1SS 


NUE 


ANUE 


0.011 


1 


.21 


Z2SS 


— > 


Z1SS 


CH 


CB 


0.002 





.25 


Z2SS 


— > 


Z1SS 


ST 


SB 


0.022 


2 


.48 


Z2SS 


— > 


Z1SS 


NUM 


ANUM 


0.011 


1. 


.21 


Z2SS 


— > 


Z1SS 


NUT 


ANUT 


0.011 


1. 


.21 


Z2SS 


— > 


Z1SS 


GAMMA 




0.785 


88 


.79 



MASS of Z2SS = 49.7181472778320 GeV 

Total WIDTH of Z2SS = 8 . 838237113852746E-010 GeV 
Lifetime = 2 . 2289512E-05 cm 



Z3SS 


— > 


Z1SS 


UP 


UB 


7815. 


.816 


4. 


.86 


Z3SS 


— > 


Z1SS 


DN 


DB 


10375. 


.482 


6. 


.45 


Z3SS 


— > 


Z1SS 


NUE 


ANUE 


6575. 


981 


4. 


.09 


Z3SS 


— > 


Z1SS 


E- 


E+ 


3426. 


.354 


2. 


.13 


Z3SS 


— > 


Z1SS 


CH 


CB 


7759. 


.725 


4. 


.82 


Z3SS 


— > 


Z1SS 


ST 


SB 


10372. 


.120 


6. 


.44 


Z3SS 


— > 


Z1SS 


NUM 


ANUM 


6575. 


981 


4. 


.09 


Z3SS 


— > 


Z1SS 


MU- 


MU+ 


3426. 


.271 


2. 


.13 


Z3SS 


— > 


Z1SS 


BT 


BB 


9638. 


.576 


5. 


.99 


Z3SS 


— > 


Z1SS 


NUT 


ANUT 


6575. 


.981 


4. 


.09 


Z3SS 


— > 


Z1SS 


TAU- 


TAU+ 


3393. 


.080 


2. 


.11 


Z3SS 


— > 


Z1SS 


GAMMA 




16. 


.500 


0. 


.01 


Z3SS 


— > 


Z2SS 


UP 


UB 


9585. 


.074 


5. 


95 


Z3SS 


— > 


Z2SS 


DN 


DB 


11596. 


.827 


7. 


.20 


Z3SS 


— > 


Z2SS 


NUE 


ANUE 


5487. 


.067 


3. 


41 


Z3SS 


— > 


Z2SS 


E- 


E+ 


4263. 


929 


2. 


65 


Z3SS 


— > 


Z2SS 


CH 


CB 


9491. 


563 


5. 


.90 


Z3SS 


— > 


Z2SS 


ST 


SB 


11592. 


.267 


7. 


20 


Z3SS 


— > 


Z2SS 


NUM 


ANUM 


5487. 


.067 


3. 


.41 


Z3SS 


— > 


Z2SS 


MU- 


MU+ 


4263. 


.800 


2. 


.65 


Z3SS 


— > 


Z2SS 


BT 


BB 


10586. 


.754 


6. 


.58 
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Z3SS 


— > 


Z2SS 


NUT 


ANUT 


5487 


.067 


3. 


.41 


Z3SS 


— > 


Z2SS 


TAU- 


TAU+ 


4210. 


.220 


2. 


.62 


Z3SS 


— > 


Z2SS 


GAMMA 




958 


.936 


0. 


.60 


Z3SS 


— > 


W1SS+ 


DN 


UB 


549. 


.331 


0. 


.34 


Z3SS 


— > 


W1SS+ 


E- 


ANUE 


318 


.880 


0. 


.20 


Z3SS 


— > 


W1SS+ 


ST 


CB 


515. 


.701 


0. 


.32 


Z3SS 


— > 


W1SS+ 


MU- 


ANUM 


318. 


.860 


0. 


.20 


Z3SS 


— > 


W1SS+ 


TAU- 


ANUT 


298 


,372 


0. 


.19 



MASS of Z3SS = 98.9230117797852 GeV 

Total WIDTH of Z3SS = 1 . 609635813320035E-004 GeV 
Lifetime = 1 . 2238793E-10 cm 



Z4SS 


— > 


Z1SS 


UP 


UB 


64341. 


.051 





.20 


Z4SS 


— > 


Z1SS 


DN 


DB 


71851. 


.289 





.23 


Z4SS 


— > 


Z1SS 


NUE 


ANUE 


799706. 


.625 


2 


.53 


Z4SS 


— > 


Z1SS 


E- 


E+ 


1257176. 


.125 


3. 


.98 


Z4SS 


— > 


Z1SS 


CH 


CB 


64144. 


.719 





.20 


Z4SS 


— > 


Z1SS 


ST 


SB 


71837. 


.734 





.23 


Z4SS 


— > 


Z1SS 


NUM 


ANUM 


799710. 


.813 


2 


.53 


Z4SS 


— > 


Z1SS 


MU- 


MU+ 


1257128 


.750 


3. 


,98 


Z4SS 


— > 


Z1SS 


BT 


BB 


70013 


.922 





.22 


Z4SS 


— > 


Z1SS 


NUT 


ANUT 


800902 


.125 


2 


,54 


Z4SS 


— > 


Z1SS 


TAU- 


TAU+ 


1239106. 


.500 


3. 


.93 


Z4SS 


— > 


Z1SS 


HOL 




11142074. 


.000 


35. 


.31 


Z4SS 


— > 


Z2SS 


UP 


UB 


70878 . 


.781 





.22 


Z4SS 


— > 


Z2SS 


DN 


DB 


113081. 


.391 


0. 


.36 


Z4SS 


— > 


Z2SS 


NUE 


ANUE 


47392 


.527 





.15 


Z4SS 


— > 


Z2SS 


E- 


E+ 


651443. 


.563 


2 


.06 


Z4SS 


— > 


Z2SS 


CH 


CB 


70611. 


.727 





.22 


Z4SS 


— > 


Z2SS 


ST 


SB 


113055 


.344 


0. 


.36 


Z4SS 


— > 


Z2SS 


NUM 


ANUM 


47392. 


.547 





.15 


Z4SS 


— > 


Z2SS 


MU- 


MU+ 


651432 


.125 


2 


.06 


Z4SS 


— > 


Z2SS 


BT 


BB 


109480. 


.492 


0. 


.35 


Z4SS 


— > 


Z2SS 


NUT 


ANUT 


47398 


.039 





.15 


Z4SS 


— > 


Z2SS 


TAU- 


TAU+ 


647318. 


.063 


2. 


.05 


Z4SS 


— > 


Z2SS 


HOL 




1357764. 


.000 


4. 


.30 


Z4SS 


— > 


Z3SS 


NUE 


ANUE 


500639. 


.094 


1. 


.59 


Z4SS 


— > 


Z3SS 


E- 


E+ 


8530 


.695 





.03 


Z4SS 


— > 


Z3SS 


NUM 


ANUM 


500641. 


.719 


1. 


.59 


Z4SS 


— > 


Z3SS 


MU- 


MU+ 


8529. 


.966 





.03 


Z4SS 


— > 


Z3SS 


NUT 


ANUT 


501390. 


.094 


1. 


.59 


Z4SS 


— > 


Z3SS 


TAU- 


TAU+ 


8321 


.920 





.03 


Z4SS 


— > 


W1SS+ 


DN 


UB 


26603 


.908 





.08 


Z4SS 


— > 


W1SS+ 


E- 


ANUE 


2805715. 


.750 


8 


.89 


Z4SS 


— > 


W1SS+ 


ST 


CB 


26414. 


.109 





.08 


Z4SS 


— > 


W1SS+ 


MU- 


ANUM 


2805695. 


.500 


8 


.89 


Z4SS 


— > 


W1SS+ 


TAU- 


ANUT 


2799321. 


.750 


8 


.87 



MASS of Z4SS = 124.350257873535 GeV 

Total WIDTH of Z4SS = 3 . 155728429046366E-002 GeV 
Lifetime = 6 . 2426160E-13 cm 



W1SS+ 


— > 


Z1SS 


UP 


DB 


11165. 


.424 


20. 


.37 


W1SS+ 


— > 


Z1SS 


NUE 


E+ 


3578. 


.382 


6. 


.53 


W1SS+ 


— > 


Z1SS 


CH 


SB 


11047. 


.252 


20. 


15 


W1SS+ 


— > 


Z1SS 


NUM 


MU+ 


3578. 


.319 


6. 


.53 


W1SS+ 


— > 


Z1SS 


NUT 


TAU+ 


3548. 


.031 


6. 


.47 


W1SS+ 


— > 


Z2SS 


UP 


DB 


7352. 


.349 


13. 


.41 
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W1SS+ — > Z2SS NUE 

W1SS+ — > Z2SS CH 

W1SS+ — > Z2SS NUM 

W1SS+ — > Z2SS NUT 

MASS of W1SS+ 

Total WIDTH of W1SS+ = 

Lifetime = 



W2SS+ 


— > 


Z1SS 


UP 


W2SS+ 


— > 


Z1SS 


NUE 


W2SS+ 


— > 


Z1SS 


CH 


W2SS+ 


— > 


Z1SS 


NUM 


W2SS+ 


— > 


Z1SS 


NUT 


W2SS+ 


— > 


Z2SS 


UP 


W2SS+ 


— > 


Z2SS 


NUE 


W2SS+ 


— > 


Z2SS 


CH 


W2SS+ 


— > 


Z2SS 


NUM 


W2SS+ 


— > 


Z2SS 


NUT 


W2SS+ 


— > 


Z3SS 


UP 


W2SS+ 


— > 


Z3SS 


NUE 


W2SS+ 


— > 


Z3SS 


CH 


W2SS+ 


— > 


Z3SS 


NUM 


W2SS+ 


— > 


Z3SS 


NUT 


W2SS+ 


— > 


W1SS+ 


UP 


W2SS+ 


— > 


W1SS+ 


DN 


W2SS+ 


— > 


W1SS+ 


NUE 


W2SS+ 


— > 


W1SS+ 


E- 


W2SS+ 


— > 


W1SS+ 


CH 


W2SS+ 


— > 


W1SS+ 


ST 


W2SS+ 


— > 


W1SS+ 


NUM 


W2SS+ 


— > 


W1SS+ 


MU- 


W2SS+ 


— > 


W1SS+ 


BT 


W2SS+ 


— > 


W1SS+ 


NUT 


W2SS+ 


— > 


W1SS+ 


TAU- 



MASS of W2SS+ 

Total WIDTH of W2SS+ = 

Lifetime = 

MASS of UPL 

MASS of DNL 

NUEL — > Z1SS NUE 

NUEL — > Z2SS NUE 

NUEL — > Z3SS NUE 

NUEL — > W1SS+ E- 

MASS of NUEL 

Total WIDTH of NUEL = 

Lifetime = 



EL- 


— > 


Z1SS 


E- 


EL- 


— > 


Z2SS 


E- 


EL- 


— > 


Z3SS 


E- 


EL- 


— > 


Z4SS 


E- 


EL- 


— > 


W1SS- 


NUE 


EL- 


— > 


W2SS- 


NUE 



E+ 2441.424 4.45 

SB 7255.331 13.23 

MU+ 2441.417 4.45 

TAU+ 2417.434 4.41 

82.2054044085650 GeV 
5.482536228723776E-005 GeV 
3.5932274E-10 cm 



DB 


389890. 


.656 


7. 


.89 


E+ 


376353. 


.063 


7. 


,62 


SB 


389994. 


.125 


7. 


.89 


MU+ 


376345. 


.625 


7. 


.62 


TAU+ 


365778. 


.000 


7. 


.40 


DB 


626844. 


.563 


12 


.69 


E+ 


381795. 


.750 


7. 


.73 


SB 


625884. 


.000 


12 


.67 


MU+ 


381794. 


.438 


7. 


.73 


TAU+ 


381132. 


.250 


7. 


.71 


DB 


1306. 


,987 





.03 


E+ 


154496. 


.859 


3. 


.13 


SB 


1267 


.973 





.03 


MU+ 


154476 


.875 


3 


.13 


TAU+ 


141892. 


.547 


2. 


.87 


UB 


2861. 


.109 





.06 


DB 


5163 


.812 





. 10 


ANUE 


41784 


.742 





.85 


E+ 


15585. 


.073 





.32 


CB 


2810 


.073 





.06 


SB 


5158. 


.582 





. 10 


ANUM 


41784. 


.945 





.85 


MU+ 


15583. 


.285 





.32 


BB 


4413 


.938 





.09 


ANUT 


41842 


.730 





.85 


TAU+ 


14907 


.641 





.30 



122.500414227970 GeV 
4.941149627879395E-003 GeV 
3.9869262E-12 cm 

259.504386750604 GeV 

263.370394294766 GeV 

17929614.000 11.16 

88574.359 0.06 

11599536.000 7.22 

131070248.000 81.57 

116.793041222755 GeV 
0.160687973351000 GeV 
1.2259785E-13 cm 

147906416.000 49.85 

140011936.000 47.19 

4365978.000 1.47 

87047.195 0.03 

1985265.250 0.67 

2332605.750 0.79 
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MASS of EL- 

Total WIDTH of EL- 

Lifetime = 

MASS of CHL 

MASS of STL 

NUML — > Z1SS NUM 

NUML — > Z2SS NUM 

NUML — > Z3SS NUM 

NUML — > W1SS+ MU- 

MASS of NUML 

Total WIDTH of NUML = 



Lifetime 






MUL- — > 


Z1SS 


MU- 


MUL- — > 


Z2SS 


MU- 


MUL- — > 


Z3SS 


MU- 


MUL- — > 


Z4SS 


MU- 


MUL- — > 


W1SS- 


NUM 


MUL- — > 


W2SS- 


NUM 


MASS of MUL- 






Total WIDTH 


of MUL- 




Lifetime 






MASS of BT1 






MASS of TP1 






NUTL — > 


Z1SS 


NUT 


NUTL — > 


Z2SS 


NUT 


NUTL — > 


Z3SS 


NUT 


NUTL — > 


W1SS+ 


TAU- 



MASS of NUTL 

Total WIDTH of NUTL = 



Lifetime 






TAU1- — > 


Z1SS 


TAU- 


TAU1- — > 


Z2SS 


TAU- 


TAU1- — > 


Z3SS 


TAU- 


TAU1- — > 


W1SS- 


NUT 


TAU1- — > 


W2SS- 


NUT 


MASS of TAU1- 




Total WIDTH 


of TAU1- = 


Lifetime 






MASS of UPR 






MASS of DNR 






ER- — > 


Z1SS 


E- 


ER- — > 


Z2SS 


E- 


ER- — > 


Z3SS 


E- 



125.486749674947 GeV 
0.296689237161868 GeV 
6.6399441E-14 cm 

259.504386750604 GeV 

263.370394294766 GeV 

17929614.000 11.16 

88574.359 0.06 

11599536.000 7.22 

131069416.000 81.57 

116.793041222755 GeV 
0.160687141212113 GeV 
1.2259848E-13 cm 

147906352.000 49.85 

140012720.000 47.19 

4365900.500 1.47 

86681.125 0.03 

1985265.250 0.67 

2332605.750 0.79 

125.486749674947 GeV 
0.296689524561566 GeV 
6.6399380E-14 cm 

263.417851693807 GeV 

312.439636958576 GeV 

17929614.000 11.17 

88574.359 0.06 

11599536.000 7.23 

130831880.000 81.54 

116.793041222755 GeV 
0.160449605025549 GeV 
1.2277999E-13 cm 

147912640.000 49.83 

140262848.000 47.25 

4347101.000 1.46 

1986060.125 0.67 

2351768.000 0.79 

125.499288741710 GeV 
0.296860418547054 GeV 
6.6361155E-14 cm 

253.087070455429 GeV 

253.391977765737 GeV 

221465200.000 67.88 
103275600.000 31.65 
1522374.125 0.47 
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109.855352062128 
0.326263181300843 
6.0380703E-14 cm 



GeV 
GeV 



MASS of 


CHR 




253.087070455429 


GeV 






MASS of 


STR 




253.391977765737 


GeV 






MUR- 


-> Z1SS 


MU- 


221465760.000 


67 


.88 


MUR- 


— > Z2SS 


MU- 


103275296.000 


31 


.65 


MUR- 


-> Z3SS 


MU- 


1522204.250 


0. 


.47 


MASS of 


MUR- 


_ 


109.855352062128 


GeV 






Total WIDTH of MUR- 




0.326263714477010 


GeV 






Lif etime 




6.0380601E-14 cm 








MASS of 


TP2 




307.130371718121 


GeV 






MASS of 


BT2 




253.439037459913 


GeV 






TAU2- - 


— > Z1SS 


TAU- 


221675248.000 


67 


.89 


TAU2- - 


-> Z2SS 


TAU- 


103215376.000 


31 


.61 


TAU2- - 


— > Z3SS 


TAU- 


1477561.500 


0. 


.45 


TAU2- - 


— > W1SS- 


NUT 




131231.641 


0. 


.04 


MASS of 


TAU2- 




109.869716964501 


GeV 






Total WIDTH of TAU2- = 


0.326499416297990 


GeV 






Lifetime 




6.0337016E-14 cm 








MASS of 


GLSS 




148.387594191548 


GeV 






MASS of 


H0L 




56.3720788464574 


GeV 






MASS of 


H0H 




313.892691024644 


GeV 






MASS of 


AO 




300.000000000000 


GeV 






MASS of 


H+ 




310.070287501751 


GeV 






MASS of 


TP 




174.000000000000 


GeV 







Cross section of W1SS+ W1SS- = 0.89911E+00 pb 



G.2 The susygen.log file 

Event listing (summary) 



I 


particle/ jet 


KS 


KF 


orig 


P- 


_x 


P- 


-y 


P- 


_z 


E 


m 


1 


!E— ! 


21 


11 





0. 


.000 





,000 


91, 


.500 


91.500 


0. 


.000 


2 


!E-+! 


21 


-11 





0. 


.000 


0. 


.000 


-91, 


,500 


91.500 


0, 


.000 


3 


GAMMA 


1 


22 





0. 


.000 


0. 


.000 


-0, 


,964 


0.964 


0, 


.000 


4 


(SusyProdO) 


11 


79 





0. 


.000 


0. 


.000 


0, 


.964 


182.036 


182, 


.033 


5 


(W1SS++) 


11 


75 


4 


-17. 


.668 


-31. 


,556 


15. 


,259 


91.096 


82 


,205 


6 


(W1SS-) 


11 


77 


4 


17 


,668 


31. 


,556 


-14. 


.295 


90 . 940 


82. 


.205 


7 


NUT 


1 


16 


5 


5. 


.178 


-0. 


,308 


-7. 


,449 


9.077 


0, 


,000 
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8 


(TAU-+) 


11 


-15 


5 


1. 


541 


4. 


,919 


14. 


,326 


15. 


328 


1, 


,777 


9 


(Z2SS0) 


11 


72 


5 


-24. 


.386 


-36 


.167 


8, 


.383 


66. 


691 


49, 


.718 


10 


NUT~ 


1 


-16 


6 


16 


947 


16, 


,397 


1 


.743 


23. 


.645 


0. 


.000 


11 


(TAU—) 


11 


15 


6 


10. 


.948 


-2 


.815 


-2. 


.115 


11. 


637 


1, 


.777 


12 


(Z2SS0) 


11 


72 


6 


-10. 


.227 


17, 


,974 


-13. 


.923 


55. 


657 


49. 


.718 


13 


GAMMA 


1 


22 


9 


-3. 


.387 


-8 


.227 


3. 


.401 


9. 


,525 


0. 


.000 


14 


Z1SS0 


1 


71 


9 


-20. 


.999 


-27 


.940 


4. 


.981 


57. 


166 


44, 


.961 


15 


GAMMA 


1 


22 


12 


-4. 


.034 


-1 


.489 


-2, 


.487 


4. 


967 


0. 


.000 


16 


Z1SS0 


1 


71 


12 


-6. 


193 


19. 


.463 


-11. 


.435 


50. 


690 


44. 


,961 


17 


NUT~ 


1 


-16 


8 


-0 


015 


1. 


.559 


3. 


.381 


3. 


,723 


0. 


,000 


18 


pi+ 


1 


211 


8 


1. 


. 160 


1 


,789 


5. 


.627 


6. 


019 


0. 


,140 


19 


(eta) 


11 


221 


8 


0. 


.396 


1. 


.571 


5, 


.318 


5. 


,587 


0. 


.548 


20 


NUT 


1 


16 


11 


3. 


.007 


-0. 


.958 


0. 


.086 


3. 


,157 





.000 


21 


(rho-) 


11 


-213 


11 


7. 


.941 


-1. 


.857 


-2. 


.202 


8. 


480 


0. 


.741 


22 


GAMMA 


1 


22 


19 





.076 


1 


,008 


3. 


.876 


4. 


,006 


0. 


,000 


23 


GAMMA 


1 


22 


19 


0. 


.319 


0. 


,563 


1. 


.442 


1. 


581 


0. 


.000 


24 


pi- 


1 


-211 


21 


0. 


,341 


-0 


,086 


-0. 


.036 


0. 
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,140 


25 


(piO) 


11 


111 


21 


7. 


.600 


-1 


.771 


-2. 


.165 


8, 


099 


0. 


.135 


26 


GAMMA 


1 


22 


25 


0. 


.863 


-0 


.245 


-0. 


.241 


0, 


929 





.000 


27 


GAMMA 


1 


22 


25 


6. 


.736 


-1 


.526 


-1. 


.924 


7. 


,170 


0. 


.000 






sum: 


0.00 







,000 


0. 


,000 


0. 


.000 


183. 


000 


183. 


.000 



Susygen has now finished . . . Bye 
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